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ABSTRACT 


A *0¢a(p5 2D) experiment was performed at the 
PIGhUME sracirrcy sine a 1200 MeV polarized proton beam. 
Data were acquired at eleven independent pairs of angles 
for separation energies from 0 to 100 MeV over a wide range 
of recoil momenta. For eight angle pairs, cross sections 
and analysing powers were extracted from the data as a 
Pureecion Of = they energy “difference: of the final state 
protons for knockout from the valence states Clidig 9's 2850/5 > 
and lds »5)- Cross section and analysing power information 
Was" ODtained= for 3°-other “angle pairs for only the 2817/2 
Siac Cv. 

The results were compared with Distorted Wave Impulse 
Approximation (DWIA) calculations which incorporated the 
Dovuncgmstateuwave functions of Elton ‘and Swift (EL-67), an 
optical model potential, wlth =spin-orbit terms, from 
Schwandt (SC-80), and a half-off-shell prescription for the 
free p-p scattering matrix element. The %£=2 states show 
strong j-dependence in the analysing power at most measured 
angles. Such dependence was predicted by Jacob and Maris 
(JA-73) and has been observed for 1 states in /!% CKi= 
S05 K 1-76) < DWIA calculations agree well with the present 
analysing powers for equal and near equal angles in the 
forward direction. The spectroscopic factors resulting 


from normalizing the calculations to the data are only ~502 
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fo ~55Z%Z of the shell model values of (2J+1)- This is 
significantly lower than spectroscopic factors resulting 
from other reactions. Generally, the kinematic regions of 
better agreement between the analysing power data and 
calculations are similar to Sines of the earlier 
169(3,2p) experiment. fies. aualys ine. power for the 
28)/,2 State near %=0 dominance is reasonably predicted 
whether spin-orbit dependence in the optical model is 
included or not. 

Ape varioudgski nematic. comditions, peaks. in the cross 
section suggest the possible presence of l1pj,;) and 1p3,/92 
Oca Pin elemeparitcleumatates at wo 21, MeV sand .~30 MeV 


separation energy respectively. 
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reaction illustrating the labelling of the 
various kinematic quantities. 


Experimental configuration of polarimeter, 
scattering chamber, and detectors. 


Logic diagram of polarimeter electronics. 
Logic symbols are defined in Figure ive 


Diagram of Nal photomultiplier tube base. 


Logic diagram of electronics located near 
scattering chamber. 


First of two logic diagrams of electronics 
located in remote counting room. 


Second of two logic diagrams of electronics 
located in remote counting room. 


Timing spectrum for start=right plastic, 
stop=every second rf pulse. Gounts are LF- 
RB (25°-67°) events with “Ca otarget ‘and 
with no decelerated beam present. 


Timing spectrum for ‘start=right plastic, 
stop=every second rf pulse. Experimental 
conditions are as for Figure III-4 along 
with the presence of decelerated beam. 


Timing spectrum for start=left, stop=right 


plastic. Counts are LF-RF (30°-30°) events 
with a +a target. 
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Counts per unit area for each of the 4 Nal 
detectors: as a function of» distance» (in 
wiresijawire spacing! c62mm) ‘from ‘the center 
of the detector. Counts are LF-RB (30°-54°) 
and LB-RF (49°-35°) events with a 40Cq 
target. Events falling outside a radius of 
26 wires (indicated by arrows) were rejected 
in the analysis. 


Distribution of events (expressed as_ the 
percentage of total) passing through y plane 
(ikea tee2o)) was eamrunetion of span. Events 
with spans greater than 4 (indicated by the 
arrow) were considered multiple hits. 


Scatter? plot (logaritthmie density’ scale)n for 
Tip versus Tap: Counts are LF-RB (30°-54°) 


events with a *%a target. Thee Gp, 2p)? Loci 
as well as the free p-p peak and reaction 
tails “are labelled. the veaction tail ‘on 
the Low -enerey side "(LB) is only barely 
visible. 


Scatter plot (logarithmic density scale) for 
Kx @@yergus = Ay/2-5 Cproportional ~*to” ‘opening 
angles in the horizontal and vertical planes 
respectively). Counts are LF-RB (30°-54°) 
events with a ‘+a tareet. The free p-p 
events (selected by software windows) and 
Chrererewo nemaining= (Dp, 2p) events “are 
labelled. Events Pala aa Watt unt the 
peetaneular region ~were,” rejected in the 
analysis. 


Scatter plot (logarithmic density scale) for 
ThoF versus Tree Conditions are the same as 


Gom Eigune IV=3. but with .all vents-pfialling 
within the rectangular window shown in 
Figure IV-4 rejected. Although >96% of the 
free p-p events have been eliminated, the 
p-p peak is still visible because of the 
logarithmic scale. 


Nal pedestal pulse height centroids for each 
detector as a function of run number. Data 


shown are for PICA uaT US: o Wi Le yi gaps 


correspond. to acquisition of CH, or CD» 


data. 
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Nal LED pulse height (minus pedestal) 
centroids for each detector as a function of 
run number. Conditions as for Figure IV-6. 


Free p-p Nal pulse height (minus pedestal) 
centroids for each detector as a function of 
run number. Data are from scattering from 
hydrogen contamination in +%Ca target. For 
each set of data, the detector angle and the 
difference between the measured pulse 
heights and the plotted values are given. 


Kinematically determined energy at Nal 

minus the energy determined from the Nal 

free; tprp ipudkserpihetightis: 10, 2° Tee OT 

Sache of the: 4) idetectors as® a function of 

Tha . Pulse heights are from CH, as well as 
in 

*0Ca runs. 


Same as in Figure IV-9 after nonlinearity 
correction has been applied. 


Scatter plot (logarithmic density scale) for 
Tip versus Trp: Counts are LB-RB (44°-44°) 
events with a CH» target. The free Pap peak 
and reaction tails as well as the L2e6(p, 2p) 


locus are labelled. 


Relative plastic pulse heights .for each 
detector as a function of vertical distance 
from the center of the event-defining 
window. Data are “from free p-p’ scattering 
from “the” CH “target. The’ "posttion or” the 
light guide and the center of the window are 
indicated. 


Counts as a function of the mass parameter 
for the RB detector. Counts are LF-RB (30°- 


54°) events with a [ae targets©* Proton and 
deuteron peaks are labelled. 


Cato 2D) cross sections as a function of 
ESUMe = (in) channels) for LF-RB (30°-54°), 
EDIF=20 MeV. The positions of the +*0Ca 
valence states (ld3,;9, 28 j;/2, 145,/2) are 
indicated along with those of possible mOCa 
lp states. The wide and narrow-fit windows 
are indicated by horizontal single and 
double headed arrows respectively. 
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Distribution of RAT's, relative numbers of 
169 and» *°¢a atoms, produced from first pass 
fitting of data. TVeweceneroia Of Che 
distribution is indicated by an arrow. Pie? 


Results of fitting LF-RB (30°-54°) EDIF=20 
MeV data. Data points and their errors are 
indicated by asterisks and bars and _ the 
overall fit by a dashed line. The ‘Ca 
(143/25 2817/2, 1ld5/2) and 1°O (1piy9s 1P372) 
peaks are shown as solid lines and dashed 
lines with circles respectively while the 
background peak is denoted by dashed lines 
with pluses (t+). The ‘energy calibration 
shift has been applied to the ESUM scale. 116 


Distribution of percentage points from the 
second pass fitting of the data. ems) 


Distribution of x? per degree of freedom 
from the second pass fitting of the data. Lis 


Cross sections and analysing powers for the 
ld3,2, lds 2, and 2s,,9 states. Results for 
the Id2 75 (circles) and Ide) (crosses) are 
presented as a function of energy sharing 
(EDIF) for the angle’ pairs indicated while 
the, 28,75 results are presented as a 
function of the free p-p center of mass 
eeattering. angle (as. calculated using .a 
half-off-shell prescription). The DWIA 
calculations for the 1ld3,5, lds 2, and 2814/2 
states are given as solid, dashed, and solid 
lines respectively while the dashed line in 
the plot of 28/9 analysing’ powers 
Gorresponds to the free p-p analysing 
powers. Liven WEA calculations “were 
performed using spectroscopic factors one- 
half the simple shell model values of 2J+1l. 122-124 


As in figure VI-l except that the spin-orbit 
term has been turned off in the DWIA 
calculations. dL 


Asein? tigure Viel @except that 2 nonlocality 
correction (method i) has been incorporated 
in the DWIA calculations. 135 


As in figure VI-1l except that a nonlocality 


correction (method ii) has been incorporated 
in the DWIA calculations. rows 
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VI-5 Plot of £ +P ie ¢ (ld 3/2) (circles) and 
—( &+1) °P AAO) Cerosses)= as a function 
of EDIF for the angle pairs indicated. 141 
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CHAPTER I 


INTRODUCTION 


Quasi-elastic (p,2p) scattering occurs when the 
incident proton undergoes a single violent interaction with 
rhremuuctear protons, thie collision “is .weuallly treated as 
free proton-proton scattering. Thus, although there has 
been some theoretical investigation of two step processes 
(KU-79), quasi-elastic scattering is usually viewed as a 
single step process as represented schematically in Figure 
fi) Currently, one of the most immediate aspects of 
(p.2p) reactions is that they ellow a check of the validity 
of the wave functions describing the propagation of protons 
through the nucleus. Alcousein, /4dcd1 ction, to ,~orfLering an 
opportunity to examine the proton-~proton interaction in. the 
Bawironnent sof era nwuctear potential, a study of the 
resulting hole states in the residual nucleus tests nuclear 
models by providing information on the overlap between the 
wave functions of the target ground state and ee ed 
states of the residual nucleus. 

A classical view of the (p,2p) reaction is as follows. 
As illustrated in Figure I-2, an incoming proton (mass M) 
which is incident upon a target nucleus (mass Mag) strikes 
a bound proton and the two protons scatter to detectors 
without undergoing any other violent interactions. In the 
shell model, each bound state nucleon is characterized by a 


well-defined orbital angular momentum and binding energy. 
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MLCU Ria 
Schematic represent atione. ot aya (fp, ec 
reaction as a single step process. 


EeUGU RE ie 2 
Diagrammatic prepresentatiio pas ote samrp fm 
reaction illustrating the labelling of the 
various kinematic quantities. 
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Thust@thes interaction fésults” in ‘the ‘knocking out of a 


moving nuclear proton which “almost” behaves like a free 


nucleon. Such "quasi-free" collisions result in a strong 
angular correlation between the outgoing protons as well as 
peaks in the summed energy spectrum corresponding to the 
binding energies of the various ings particle states. 

A characteristic of quasi-free scattering is that the 
tareet “nucleon tan “haverYan Yetfective initial polarization 
oaeere eaction In asymmetric kinematic conditions, the 
reaction may be localized. For example, below 400 MeV, the 
mean free path of protons travelling through nuclear matter 
decreases strongly as the energy decreases. Consequently, 
events in which one proton emerges with a large energy 
compared to the other may tend to Orieinate on“the side +of 
the nucleus in the direction of the momentum of the lower 
energy proton. This results ‘ina’ ‘spatial “localization of 
the reaction. @Thistis *distince irom the®localizatton of 
the (p,2p) reaction due to radial peaking of the bound 
state wave function. Surface localization does not occur 
Poreecie 23.775. Suate™ “wave furction “whirene@"is mainly 
concentrated near r=0.) A careful selection of the angles 
and energies of the outgoing protons determines the 
magnitude and direction of the momentum off sthe* struck 
particle which can couple to a localized enhancement of the 
reaction and this leads to a well-defined sense of orbital 
angular momentum. Through nuclear spin-orbit coupling, the 


struck proton is effectively polarized in a direction 
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normal to the scattering plane. in particular, protons in 
States of the same & but different j will be polarized in 
opposite directions. This, coupled with the fact that the 
free p-p cross section for the scattering of protons with 
nanallel. spins, is slarger thah that for antiparallel spins, 
leads to a. j-dependence in the (p,2p) analysing power. 
This dependence was predicted by Jacob and Maris (1A ~=73.) 
and first observed for the %=1 states of 169 (KI-80,KI-76).- 
Therefore, CPp..2 Pp.) could conceivably be used as a 
spectroscopic tool to distinguish states of the same & but 
ditferent. j. 

AW auantitative, picture of Gs 2p) Lineludes. the 
effect of nuclear distortion of the proton wave functions 
and is formulated using the Distorted Wave Impulse 
Approximation (DWIA). In Appendix I, the DWIA description 
is discussed and details are given about the manner in 
which the j-dependence in the analysing power arises. In 
the present ;context, the struck ‘particle polarization and 
hence the analysing power are more promising than the cross 
section in testing DWIA. This. arises . because .of_.the 
particular insensitivity of analysing power to the bound 
state wavefunction and the spectroscopic factor. Also, the 
analysing power is less sensitive to the optical potential 
than the cross section. 

From the earliest experiments to the present, analysis 
of (p,2p) experiments has consisted of organizing the data 


either on the basis of the summed energy of the outgoing 
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protons, on the recoil momentum of the residual nucleus, on 
the outgoing angles, or on some combination of these (KU- 
PET TA=G69UPEY766 , WE=557CH-52)% Eanly interest. in Cp,2p) 
reactions focussed on attempts to study the momentum 
distribution of ound *nueleons by® measuring” "energy 
distributions for one outgoing ar ceed over a variety of 
scattering angles (CL-52). In the first (p,2p) coincidence 
experiment’ *(CH=52);,°tangul’ar’® correlations *for SLi were 
measured about an opening angle 90°. The degree to which 
bound nucleons exhibited free nucleon (quasi-free) behavior 
was examined in these _ studies. ave re coaneraence 
experiments which measured the energy of one E€WE=55 ) © or 
both (TY-58) outgoing protons gave some of the first direct 
evidence of the existence of single particle shell model 
states. 

Several features of quasi-free scattering contribute 
to the Sinterese *in® this -reaction. The kinematics of the 
(p,2p) reaction allow the recoil momentum to vary over a 
wide and useful range from zero to several hundred MeV/c. 
Consequently, it should be possible to distinguish between 
states of differing® %,° particularly 240 ° states (whose 
momentum distribution, and hence cross section, is at a 
minimum at zero recoil momentum - as explained in Appendix 
I) from &=0 states (whose momentum distribution is at a 
maximum there). In addition, the same binding energy and 
recoil momentum region can be examined with a number of 


different kinematic configurations. This redundancy allows 
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one to examine the interaction process for a given nuclear 
state under a variety of conditions. Comparisons with a 
DWIA formalism can then be made. 

During the quasi-elastic reaction, the binding energy 
of the struck nucleon together with the recoil momentum of 
the residual nucleus makes the 2 body t=natrix, half-ort- 
shell (RE-70). In addition, the distortion suffered by the 
incoming and outgoing protons causes the t-matrix to be 
fully off-shell (RE-/3). However, if this latter effect 
can be ignored and the DWIA is a sufficiently accurate 
description of the (p,2p) reaction, it may be possible to 
extract oOff—shell top Matrix information which ‘can 
then be compared with the predictions of various phase- 
equivalent (phenomenological) potentials Ca 1-767)". Such 
potentials are equivalent on the energy shell and _ the 
hope is that one can distinguish among them by calculations 
of off-shell processes. Unfortunately, to extract nucleon= 
nucleon information one must perform accurate calculations 
of the effects of distortion, which can be quite large. 

Tt *has been suggested that certain experimental 
geometries might be Dartutcurarly propitious “for the 
examination of off-shell effects (I10-78,DE-/9). With such 
geometries an attempt would be made to utilize the variety 
of kinematic conditions possible for the reaction. The 
expectation was that the sensitivity of calculations to 
distortion and the bound state wave functions could be 


reduced by holding the outgoing proton energies constant 
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while the initial state of the nucleus being examined was 
fixed aby mestricting) ,thepyvaloe.of;)the recoil momentum. 
However, calculations comparing off-shell behavior in 
(p,2p) geometries expected to maximizes hsuch.« differences 
have shown disappointingly small sensitivity to the 
potential used (MI-/9). 

From the above discussion, the (p,2p) reaction is seen 
tObn Ot fe information on. nuclear structure, off-shell 
effects, and propagation of protons through nuclear matter. 
There are, ,.of; course,,ya number ,of, other reactions which 
give information complementary to Ehatesoptatined ping Cp, Zpde 
Through knockout or pickup reactions suc s<asybiGp, pn), 
(eker pi)necGc, OpjaweGpd)» (d, *He), andes Y¥pot inftormmationron 
nuclear structure through the removal of a single nucleon 
is acquired. The knockout reactions (p,pn) and (e,e'p) are 
fairly well described by a single step DWIA formalism at 
intermediate energies. bivceutthess (up a2p)alreaction; these 
offer considerable kinematic freedom including wide 
waar iraits. 6 nan Oteet hele te.coig, momen tum. Experimental 
investigations of (p,pn) reactions can be designed to 
acquire (p,2p) information simultaneously (JA-79). Such 
measurements are expected to yield information about bound 
neutrons similar to that from (p,2p) reactions concerning 
bound protons. 

The impulse approximation can be checked by finding 
quasi-free processes which have nearly equal distortion and 


differ primarily in the knock-out process itself. An 
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example of such a test was an analysis of 12c(p, 2p) and 
12¢(p, pn) datauGMA=4 9). By .taking the ratio /of the.cross 
sections for each process, the effects of distortion nearly 
cancel. A comparison of the ratio of the measured cross 
sections with the theoretical predictions serves as a 
check of the impulse approximation. 

Information about the momentum distribution for the 
deeply bound states of several nuclei including *lCa has 
been acquired through (e,e'p) measurements (MO-/76,NA-/76,AM-— 
66). Compared-to (p,pN). (e,¢'p), with only one strongly 
interacting particle, might .be expected to suffer less 
nuclear distortion and the interaction is presumably well 
understood. However, cross sections tend to be several 
orders of magnitude lower and thus harder to extract from 
background. Enticines evidence. yfor ts. and. 1p single 
particle state strengths in 40¢a has been published (MO-76) 
although better statistics would be desirable. The (Y,p) 
reaction can also excite hole states and offers advantages 
and disadvantages similar to those for Ce peop), but. the 
range of experimentally accessible recoil momenta is 
restricted to values >300 MeV/c (MA-75).- Apart.fcom ithe 
factiathabs datmuigeuddimitedseidue. to..the difficulty; y,of the 
experiments, questions about the reaction mechanism make an 
interpretation of the results diet cilities AA te C.e nt leo( y,p) 
experiment suggests that one-step processes are inadequate 
to account for the observed differential cross sections and 


that the inclusion of two-step processes improves agreement 
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wth “the data '€MA-7'7')'. 

(a,op) quasi-elastic scattering is strongly surface 
localized compared to (p,2p) due to the strong absorption 
of *taS particles: in nuclear matter. Some recent (4, ap) 
measurements on °Li and !%F are in fair agreement with DWIA 
calculations and, in the case of Staiy, yield a spectroscopic 
factor consistent with that derived from a DWIA analysis of 
SLi(p,2p) data” (NA=79)s Measurements “of! *(a,ap) and (p','2p) 
on +a are currently being analysed (RO-80). 

ThReeeneutronyrpLlckupy reactd’on ,PCpiid)4° provide's 
information about high momentum (>300 MeV/c) components of 
bound neutron wave functions. Recent Heo Clat(?p tds) 
measurements show evidence for a number of single particle 
states up to excitation energies of 18 MeV (LI-80,AL-79,SM- 
Toy. In “addition; this reactifon*is potentially’ useful ‘in 
probing deep hole states. With the advent of polarized 
beams, (p,d) analysing power measurements have been 
shown to distinguish strongly between states of the same & 
Hut’ tditte rent £4** (HU=s80, AL=7:9', CH-67)< High momentum 
components are also examined in the (d, *He) reaction. 
Stategurot Bditterenteres ane Sreadilyy*distinguished-* by’ the 
angular’ %distributions “of ‘their cross )*sections GKR=71 = 
High resolution *0Ca(d, 2He) measurements at 52 MeV have 
provided detailed spectroscopic information from the ground 
state to excitation energies of ~10 MeV (DO-/6). 

The current investigation of SVeaC ph S.2 Dp) iS san 


extension of the previous 169(5,2p) experiment (KI-80, 
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KI-76) and offers new possibilities. The choice of 169 and 
+ Ca targets follows from their reputations as good shell 
model nuclei which should provide relatively pure single 
particle, states. 

The separation between any two of the three valence 
states of ‘Ca (dee. west p8tadicgy) d isedessacthame3 dMeV 
making them experimentally more difficult to resolve than, 
for example, the states of 169 which are 6-5 MeV apart. 
With ‘improvingy-energy »resolution,, more) recent +*Ca(p,2p) 
experiments (TY-66,JA-69,KU-/1) have tended to emphasize 
the nuclear spectroscopy aspect. A recent +* Ca(p, 2p) 
experiment (RO-78) attained a resolution of 350 keV but 
examined only the ld3,, and 28,9 states. The best 
resolution attained in a (p,2p) measurement which examined 
all three ‘Ca valence states was 4 to 5 MeV (JA-69).- 
Tive,erca, goal pote this,+expeniment,wais to, jathtainr thes best 
possible beam and detector resolution in order to fully or 
partially!) -resolve’ all .thesvalence) states. A knowledge of 
excitation energies and relative spectroscopic factors for 
the removal of protons from a particular single particle 
state can greatly assist in the analysis and interpretation 
of nucleon knockout experiments where the resolution may be 
inadequate to resolve individual peaks. Such efactors? for 
*0Ga have been provided by DWBA analyses of (p,d), (d, 2He), 
and (dehy experiments CDO-76,MA-/Z). and were helpful in 
this investigation. 


In the present experiment, both the cross sections and 
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analysing powers for the three valence states of *0Ca were 
measured over a large range of recoil momenta. With the 
kinematics chosen so as to achieve a large effective 
polarization of the struck proton, the j-dependence of the 
analysing power of %=2 states was examined for the first 
time. The degree to which DWIA could reproduce the data 
and, “in /particular, the®*extent™ to’ which the» inclusion of 
spin-orbit distortion affected agreement of the predictions 
With the “data were -of"mayor interest. This experiment was 
the first experimental investigation of Ca r2py ontt’a 
medium mass nucleus and was also the first investigation of 
the ability of DWIA to reproduce such data. 

In previous 40Ca(p,2p) experiments, James (JA-69) and 
Kullander (KU-71) both reported evidence for lp and ls 
strength in we OF The current experiment also offers the 
prospect of examining deeply bound nuclear. states and 
provides the possibility for distinguishing between states 
such as lpj,,;5 and 1p3,; 2 on the basis 4o8- Ehe® Dehavior ‘oft 
their analysing powers. The measurement of binding 
energies and widths of deeply bound states provides 
important data for = comparison ‘with’ nwelear mat Cer 


calculations from 2-body interactions (C0O-/74,SP-72,SP-/1).- 


12 


a] r ¥ 7 : wie 
Whe = , ee 

severe sonetatv swore 4hs.tet annwes witeyiags — 
p : 


a> wt ,eseoeon Ltwees Ve -agran wogaant “so ss¥0 shysondon 


: RP I : 
4s wunwel » sevetdtod.§ 3 s oF seeds pskvedinial. 
aa 


a ae : -* : oe 
e¢ ‘ bef ofS ,fe2eTq SoMa7e CHS 70 rolsost salto: 


a 


4 +4 L @ r, ‘ “ALanu ‘’s aO'7' —j 5 9 ae sewag joleylias 


pond 92> AIW ASdte oF Yaob wit 
> I bs 
uw } 2 c Oa < #2 “twee os 744 - a2 
1 as na * ‘ane an 
“4 4y ch “zAaod & @¢ twlath théve~-ale 
: ; - ¥ 
; nd 7 7 hn 
ae i -SsevSin2e 3OSn “To Sete £9551 
: ae e 
7 ~~ £ 
ag bsanFigesvo)  Lappenlreges - 20977? 
a. ™ Ff 
5 x a 7 : ‘ 
sve! sedli 24 OCS& ZkW SAP GeOupSur 822m wii PbS or 
A” a = iy t 
Coe-% *2LSDOTOSI 7 prude? to vititve 4 
— . 
~~ 2 
‘ fa 
Tr t “pp % H rs 8 Le re ie aud} ay oj 
di + 


> vinehEes 4 fraqos tdae ~(i 4-05 ~ teboeti E 
eo te | Bay 


Tre ite tres ent -.aa* nf dsanusae 


GH I4 Se wIe0 | bogntzest ret % Dads ettg> ona etal iet. 


f _ ; — 
ja} io : ' Rea he ore Vy el — Howe 
Au 
pa 2 i J i a neg, ony i aw tee (ome ae 


7 oe a 

vsatJy pound @igast te aeehiv: Gare eelgren ; 

: 7 * a mS ; 
: Yee how 4 wk foclvcg@ad’ ae wath. 20837409 


_ ¢ 


o>yY s90)s peresal «Boek wT Antsy 
as, Vn 


CHAPTER II 


EXPERIMENTAL DETAILS 


2.1 Overview 

The experiment was performed at the TRIUMF facility 
located near the University of British Columbia. The 
TRIUMF cyclotron delivered a 200 MeV polarized proton beam 
Wipe polarization of 7-0 .o. > -0./5: and with an energy 
resolution better than 1&2. The-cyeclotron accelerates H™ 
ions and a proton beam is extracted by removing both 
electrons from the ions by means of a stripper foil. The 
proton beam has <7 ns pulses separated in time by ~43 
ns. During the experiment the average beam current was 
kepee at 0.5 to, b.0 nA. 

The experimental configuration is illustrated 
genenatctcally in Figure Lil. Four detector telescopes 
were mounted on booms which could be positioned 
independently around a circular scattering chamber with an 
angular precision of TORO. Beam polarization was 
measured by means of a polarimeter located upstream of the 
scattering chamber. Beam current was monitored both by 
means of the polarimeter and an ion chamber located 
downstream of the scattering chamber. A more detailed 
deseription of various components of the apparatus 
follows. 

Dee seolarimeter 


The polarimeter used in this experiment has_ been 
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FIGURES ers 
Experimental configuration of polarimeter, 


scattering chamber, and detectors. 
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described in detail in (GR-79,MC-78). Free p-p scattering 
Gate ae CH es target) withs an average thicknesa of 
3.83+.38 mg/cm* was used for monitoring both the beam 
polarization and beam current. The polarimeter detectors 
are illustrated in Figure I-11. Ons either. side--of -tChe 
heamlinesa “pairvor detectore at 17°. (Ll*L2 or Ril *R2) acted 
as a telescope to define the solid angle. trapvle 
coincidences (L1*L2°*L3 and R1°*R2°*R3) are called prompts 
(consisting of reals and randoms) and are counted. To 
measure the contribution of random events, coincidences 
between these same sets of detectors but with L3 and R3 
delayed by 43 nanoseconds were also counted. This delay 
ensures that the protons come from different beam bursts 
and hence different interactions. The rate of accidental 
events from one beam burst is assumed to be the same as 
that between two consecutive beam bursts. A schematic 
diagram of the polarimeter electronics is shown in Figure 
i=? is 

We denote the number of prompt coincidences L1°*L2°*L3 
(RP*R?2-R3)as POL[L] (CPOL[R}) and the number of accidental 
Potuctdences) as POL[EACC] CPOL[RACC]). The bean 
polarization ) and number of “incident protons (NE) “are 


given by 


1 POL[LRE] - POL[RRE] 


Giimwee seam Polarization == PAP © POL[LRE] + POL[RRE] 


GROWS Oy: 


Gis 2) NTs (POL[LRE] + POL[RRE]) ° eo , 
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Logic diagram of polarimeter electronics. 


Logic symbols are defined in Figure BES RU RCS ke 
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where POL[LRE] 


POL[L] - POL[LACC], 


POL[RRE ] POL[R] > — POL [RACC], 

PAP is the polarimeter analysing power which equals 

GS$30020 01S" “at” "2007 Mev, 

PC is “the cConvereion cactor from “polarimeter 

counts to nanocoulombs which equals 350425 counts 

per nanocoulomb at 200 MeV. 
Beam polarization for spin + (¥) typically ranged from 0.6/7 
PorrOee oe (-Oso! tor -O-/00e during’ ithe ‘experiment. The 
polarimeter served as the primary beam monitor in the 
experiment. 
225° Lon Chamber 

Readeetailed “deseriptiom of “the * ion *chamber, "bts 

operation, and calibration appears in (MC-78). in bDrier, 
the proton beam passing through the ion chamber ionizes 
helium gas that is continually flowing through the chamber. 
The resulting charge is collected by means of high voltage 
planes in the chamber and a digital signal proportional to 
the amount of charge is produced. During’ the expetiment 
the voltage was maintained at -300 Volts. If the number of 
ion chamber counts is represented by IC, the number of 


incident protons (NI}) is given by 


IC *6.2418x10 ? 


CLI 3} NI, ee Te Cee? 
where the conversion factor from ion chamber counts to 


nanocoulombs (ICC) is 156+t11 at 200 MeV. 


The ion chamber served as a secondary beam monitor, and its 
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consistency with the polarimeter is discussed in Chapter 
IV. 
2.4 Targets and the Scattering Chamber 

Asficitewlan scattering? chamber with 4a 7radius of°'25' ‘cm 
was used. The window of the chamber is 5 cm high, allows 
the passage of scattered particles from 122 O69 168° ton 
elitheritside off the® beamline, ‘and? «consists ‘of § 127 um? ‘thick 
Kapton. The targets used in the experiment were: CD, 
(58 mg/cm*+5Z), CH, (50 mg/cm7+5%), a? naturadevecaoltareget 
(@0mese mefenm225e2 of whith 99659778 1680 ’Caland we shail 
refer to it as the ‘Ca target pytand fabeZnsS pseintilblator 
screen for monitoring the position of the bean. They were 
mounted, one above the other, on a vertical ladder which 
can be raised or lowered remotely so as to position the 
desired target in the path of the proton beam. The target 
drive mechanism is mounted on a plate which sits on the top 
Ore tiemscatveringr chamber. lie CD’, target visecircular with 
geod bameter Sofie S few, *+while*jthe**other Yeargetsi*are™’*5 cn 
Square. 
2.5 Detector Telescopes 

The four booms surrounding the target chamber were 
designated: Lete Cron cet Le) gneki shti FrontecCRF) 5? Left >Back 
CUB), andaktent (ackt (REY, where (left oreright "refer to the 
side of the beamline looking downstream. This labelling is 
diggs trcatedtaint Pigurestli-i1s Also shown in the Figure is 
the@@arraendsenmentotott ithe’ components )“ofsethe® detector 


telescopes on each boom. Each telescope consisted of two 
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multi-wire proportional counter (MWPC) planes to provide 
position information in the horizontal (x) and vertical, (y) 
@arectianstrantbiniplastic counter, and a NaIl(T%) detector. 
The front (back) wire planes were 134.1 (115.2) cm from the 
target and were 64(x) by 128(y) (64 by 64) wires. The 
solid angles were defined using these planes. The 
construction and operation of these multiwire detectors is 
reviewed in (MC-78,CA-75). 

Each multiwire chamber consists of a plane of high 
voltage wires kept at -5000 Volts, a.sense plane of wires, 
a second high voltage plane, a second sense plane, and a 
final, high yvoltage.,piane. The sense wires are 16 um in 
diameter,. spaced 2,ammleapart;-cand’ made}.of, gold-plated 
tungsten. The high voltage wires are 102 tum in diameter, 
spaced 1 mm apart, and consist of beryllium (1-94) and 
copper. The wire planes are enclosed between sheets of 
25-42 thicks,Kkapton, and ethe.4./f cm. spacing, -between,)| the 
Kapton sheets is continuously flushed with a gas mixture 
(75% Ar bubbled through methylal at O°C, 24.5% isobutane, 
and 0.54 freon). 

AA phas Ghe giacinitil Latiom.,s 12.8 cm: squaref andr«9:635 
60742/5) 9 cayethiick awas,)positioned: behindajeach, front (back) 
MWPC. Bagheaot etheseredetiectorns ~donsists ~of«.NE1LO 
Seintillatote and isisconnected ~to,ja plastic- light guide 
which. directs. the light flashes;through a 90°, bend to a5 
cm RCA 8575 photomultiplier. The scintillators acted as 


passing counters providing timing and energy loss (dE/dx) 
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information. The energy deposited in these detectors was 
used in the calculation of a mass parameter to separate 
deuterons from protons. 

Nal detectors served as stopping counters to measure 
the kinetic energies (T) of the protons. A detailed 
description of these detectors and their efficiency appears 
im) (GCA—77))i. eines taree quantity of light. generated) by 
detected particles in intermediate energy experiments 
commonly results in a gain shift in the photomultiplier 
tube with detected flux. This problem arises because the 
current flow in the photo-cathode induces a voltage drop 
MeLaasceite which distorts ithe electric field collecting the 
photo-electrons. The electronics of the base of the 12.5 
cm photomultiplier tube were specifically designed to 
minimize this problem and are sketched in Figure II-3. The 
4 tube bases were fed through a common 10 uF capacitor kept 
at +1600 V. Normally, the outputs from the photomultiplier 
are passed through a pulse shaper before integrating with 
an analogue to digital converter (ADC). The nearly 
exponential tails of such shaped signals are quite long (10 
us) and the possibility arises that the tail of a previous 
pulse may overlap a pulse of interest. Also a second pulse 
may arrive during the integration of the first pulse. 2a 
either case the undesirable result is that the total 
integrated signal contains contributions from more than one 
pulse (a condition referred to as pile-up). To overcome 


this concern, devices called pile-up monitors which note 
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FIGURES Ll—3 
Diagram of Nal photomultiplier tube base- 
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the arrival of two signals within a specified gate (length 
of time) were used. Pile-up monitors with inputs from the 
plastic scintillators recorded instances when two signals 
frome plastic wdetector were, sufticiently close in, time 
that the corresponding Nal signals would overlap. Signals 
from the plastics were used as they are sharp Gol0 ne) and 
thus allow the pile-up monitor to resolve between signals 
close in time. Beam current was kept low enough that pile- 
ups did not exceed 20% of the signals processed. Pile-ups 
were rejected in later analysis of the data. 

In principle, the long exponential tails from the Nal 
tube base can be brought back to zero voltage (the baseline 
voltage) quickly by adding an attenuated, inverted, and 
delayed component of the original signal to the original 
signal. Such a cancellation was attempted by means of a 
clip line connected to the tube anode in parallel with the 
signal cable. Age showneimerteure LIl-3, ,the clip line -was 
Serminated swith a» .0.5. uF capacitor. and a resistance (15 
Ohms) which is lower than the characteristic impedance 
of the cable (50 Ohms). The resulting output signal was 
directed to an ADC with no pulse shaping and could be 
integrated in approximately one-tenth the time. Although 
it was expected that such clipping would decrease energy 
resolution while allowing significantly higher event rates, 
puntie Of tiie .electronics at the start of the’ experiment 
indicated that the resolution remained constant whether the 


terminators were connected or not. Subsequently it was 
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decided to proceed without the terminators and to live with 
longer integration times and necessarily lower event 
rates. 

Each Nal(T&%) crystal and tube base was encased in an 
antimagnetic shield (Mu metal) which sat in an iron pipe. 
This pipe was, in turn, mounted in a-lead pig to shield the 
Nal detector from background radiation. The Nal detectors 
were positioned 2.5 em "behind “the “plastic scintillators’. 
Broeetaetit) “crystals are cylindrical” with a -“I2.7" cm 
diameter and a 7.62 cm thickness capable of stopping 150 
MeV protons. When the front detectors (LF, RF)” were 
Operated forward of 29°, a 0.631 cm thick “copper plate was 
inserted behind each of the front plastic scintillators to 
degrade the proton energies sufficiently to stop them in 
the Nal detectors. Taped to the detection surface of each 
Nal and plastic was a light-emitting diode (LED) driven by 
a pulse generator. The pulser also fired one wire in each 
of the MWPC planes. Thes response Or the Nal detectors To 
the pulser signal had an ~50 ns rise time and a fall time 


Of =~30.0) Ds). 
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CHAPTER III 


EXPERIMENTAL PROCEDURE 


Lomepnocess (the Untormation from’ the detectors a 
combination of fast NIM and CAMAC electronics was employed. 
The electronics were divided into two sections: one near 
the scattering chamber and detectors, the other in a remote 
counting room. A Honeywell H316 computer located in the 
counting room received event information processed by the 
electronics via a CAMAC network. A. data acquisition 
program, ORION, controlled -the data processing in the 
Honeywell. Amedetaivede cescuiption of the “structure and 
eperation Of this data acquisition system has been given in 
museriLestof internal reports (RO-/4,R0-75,R0-/5a,RO-77).- 

Amecchenatic  aiacram 1tlustratinge. the electronics 
located near the scattering chamber is shown in Figure III-l 
while those located in the counting room are ‘shown in 
Figures I1l1-2 and ITII-3. A coincidence between a-left and a 
right plastic detector generated an EVENT FAST signal. This 
requested the multiwire controller to preserve any multiwire 
information in a set of memory registers (this action is 
called a strobe) and also caused the EVENT FAST logic unit 
Pee pomeracchedenso,  tiiat sno. other coincidences could’ be 
registered until a CLEAR signal was issued. The electronics 
in the counting room determined whether the event was to be 
processed further. A coincidence between a left and right 


plastic and Nal detector telescope generated an EVENT SLOW 
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FIGURE III-1 
Logic diagram of electronics located near 
scattering chamber. 
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FIGURE =bil +3 
Second of two logic diagrams of electronics 
located in remote counting room. 
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(MASTER GATE) signal. This latched the MASTER GATE and 
signaled the system to transmit the MWPC information to the 
computer through CAMAC. If a MASTER GATE was not generated 
Mithin 1|* We of the strobe, the MWPC controller would 
generate an internal CLEAR which cleared the multiwire 
memories and disabled the EVENT FAST latch so that further 
coincidences could be processed. If a MASTER GATE occurred, 
the system processed the information until it was ready for 
the next event at which time the computer issued a COMPUTER 
NOT BUSY signal. This generated a CLEAR which was sent to 
the MWPC controller and which released the latches on the 
EVENT FAST and SLOW logic. 

Following a MASTER GATE, CAMAC ADC's were turned on for 
wteixced) length of time Ca GATE) to integrate .the signals 
Heomer thew) Lasticm ana “Nal “detectors, and time-to-digital 
converters (TDC's) were started. The ADC's were multi-input 
(12-fold) charge integrating types with 256 pC corresponding 
fo “mull scale-(41024 channels). The ADC gates used were 400 
"Ws wide. The constant fraction devices shown in Figure III- 
Pepcrovidea a) Shatp timing signal when the height of the Nal 
signal had risen a pre-specified fraction of the full peak 
height (0.2 in this experiment). The widths of the pile-up 
gates were 9.0 to 9.4 us. There were separate sets of ADC's 
ander Dotseetor the Left. and right detectors. To ensure that 
aeelett, (richt). detector gated. of started the left. (right) 
ADC's or TDC*'s respectively, the: sicnals-to, initiate, these 


processes were derived by ANDing the MASTER GATE signal with 
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a suitably delayed signal from the left (right) plastic. 
This is called retiming. 

A TDC unit has’ a number of inputs each of which tis 
stopped by a separate signal and all of which have a common 
start. In addition to using timing signals from the 
plastics and Nal detectors as stops, every second signal 
from the radio-frequency (rf) oscillations of the cyclotron 
was employed as well. The resulting timing right versus rf 
Pstare=-cient , stop=every second ri pulse) spectrum 
encompassed two beam bursts. A typical spectrum is shown in 
Figure III-4. If the cyclotron becomes slightly mistuned, 
some of the accelerated H ions that pass the stripper foil 
ain Out of sphase @ with. the. rf by ~180°, experience 
deceleration, and are subsequently extracted by the stripper 
foil resulting in a decelerated beam. Data acquired during 
periods having large amounts of decelerated beam have a 
spectrum such as the one shown in Figure Pil 5 As 
discussed in Chapter IV, decelerated beam data had _ no 
objectionable features and was not treated differently from 
other data. One of the stopping signals on the left (right) 
TDC consisted of a suitably delayed signal from the right 
(left) side. A typical Jeft-right timing (start=left, 
stop=right plastic) spectrum is shown in Figure III-6.- The 
first and third peaks are random coincidence events 
consisting of protons from successive beam bursts. The 
middle peak consists of prompt (real+random) events (from 
the same beam burst). 


The DCR. units shown dn° Figures JIl-—2-vand> 21-3 are 


35 


Biteve sxeeohlonies #bhwes  o90 etasa 
By 7 = 


inal to 4 » veut Yo isedoaun > s @2eak. thaw 20f A. 


aad . a : * 


ny 

oe ie oe 

Gywi« 3 i238) a3 no077 fengze bevalsd ¢qidesios i 
: ww 


- ,getea¥es Matias af efdt 
o. a ee 


fi 7 = 
} ’ 


2 . | Aplea)4 “op sey. oa oorcintsn- al A 398 


: Pe | 5 o 142 soedab “$ ut ani #¢ 
1 ~ 1 
climifloes (23) xanechevi-pitep nas 
- a peri — ‘aly 
4 ¢ goto lweer aed ‘ hOeyol aoe 


- Sa 
a) be tug 5 bitlG tarp uteVeRRo4 , t& etre Ie78g- 


-~ ¢ = 
~ 
- 


) : : a 


ava A: foe el cad aay dreary rooay 
a ’ 7 

) a “2% oensad wereels ys “ena: SI -*=135, sougee 

ji a3e ans 2 5 ; aw Par “ti, Baterpetossa- ols Jo eros 


os r _- ; = i P - 7 = ¢ , 
7 ( *z 7 fe ~ 4 5 cat SiS is 4 See61'C +6 Fue q ite 


sae boicftu3xe6. visnsupepdue 68 bore feds sse bane: 
= a - me. 


ri aed bet, 3wiesel p ML 202088885 
j 4 a 
Esat s ei OHvuvge 2a 781 unévad 2hots “7 
— 


: ; nw ade $60 6 6SH2 (8&6 here eeoRge 
¥ : 5 


roe sS2°0Gn env “On Aato2en t. aivenp si 

% a2 ty ; 2 

%) 2tel 1 ho BELeRQe 20 1QqQ e732 Sha TO Bal 4 es To. 
: ; - — Lo 

1 lenpple bhevyatoh ~ignidoeie ve baseleags 


[esvei¢ goigss piatIn~siss fr. olfgys A athe (ts 


i dl 
-8*ITTtl g#7ugs? wh werpsta oa! vuiaeeigil ssiaasy 
— 1 


gi pier mnat ovix as yous ears a 
» : - ; 
A ; 


I 7 


: Pe oe : rand 
2 hie ai ae es ae a 


; Bo. 
oh, : 45 
7  <) 
a ' ps 
T opts 
, ® wf 
‘ “a Py : 
— SS a 
™ 
i 
®~ETS SAUDLY } 
id 2554 ° onigit 
- » = 


rc 
& 
+ >. 
e 
7 
& 
Q 
te 
os 
tu 
bow 
€ 


tt 
i] 


— at. 7. a ; Pea — ‘ae —— A 


os 


- rt err ; ' 
¢-I1T! « ee & t 
5 


y 
73627), tust725¢80 gnimti ; 
x bbnosts: Vrsve"“aois 
i wa OTe anolzgtbhnes 
dq jsoaseszg sf dite 
7 / 
io BAD 
’ ( 
ati ——_———_— = 
| a 
if) —T ro 
Le! PRUNF 1D. 
‘ ‘$ ™ % 
rr A a cir ton 
i] (RE Th. at =i j 
' rn ? 
Ht 
> i} Ne 


[ 6-JIzi szvD1% 

, ls32 to eid: ee on lm@iT -- - 

="00)-GH- 7. sts & ePry _«alszeela 

ister |)! eee 
RAND ONS 


FIGURE III-4 
Timing spectrum 
stop=every second rf pulse. 
RB §(@25 767 =)ee eventsmewlen *0Ca target 


with no decelerated beam present. 
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stop=every second rf pulse. Experimental 
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CAMAC digital coincidence registers (or bit pattern units) 
consisting of 24 bits. For each MASTER GATE, the contents 
of the register were sent to the computer with each bit 
representing the status of its input. Table I shows the DCR 
assignments used in the experiment. To record the total 
number of counts from certain devices over the course of the 
experiment, CAMAC scalers were used. The quantities scaled 
in this experiment are listed in Table II, and the positions 
Weeecome Bot these- scatters. are indicated in Figures I11I-1, 
Pee—-2, and Lii-3. 

While the computer is processing a previous event, the 
latch on) the MASTER GATE could result .in,»the loss of new 
events. This computer dead-time must be taken into account. 
In addition, the electronics preceding the MASTER GATE could 
contribute to the overall dead-time (although this effect 
is usually negligible compared to computer dead time). For 
the purposes of monitoring both the electronic and computer 
dead-times and to present a constant amplitude signal for 
calibration, pulser pulses, which were scaled (1 pulse for 
every 2°) L+R polarimeter counts, were fed to the plastic 
and Nal LED's. 

The information received by the computer was processed 
and, for each MASTER GATE, an event record was written on 
magnetic tape. A list of the contents of an event record is 
shown in Table III. In the table, LAM PATTERN. refers to. the 
CAMAC look-at-me signal issued when a CAMAC crate required 


the attention of the computer. The events written on tape 
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TABLE I 


DCR Assignments 


DCR Number Assignment 
0=5 Not used 

6 LF (Nal *P1l) 
p LB (Nal *P1) 
8 LBRF 

9 LFRB 

10 LFRF 

bik LBRB 

12 RB (Nal °P1) 
13 RF (Nal °P1) 
14 Pulser 

iS Beam on 

16 Pite Up LF 
17 Pile Up LB 
18 Pile, Up RE 
19 Pile Up RF 
20 Spin Up 

Zl, Spin Down 


22-23 Not used 
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TABLE II 


Parameters Scaled and Stored in Header and Trailer Blocks* 


Scaler Number 


OONA UN FWHM Fr 


Parameter 


Ion Chamber 

Event Slow Coincidence Presented 
Master Gate 

Elasped Time (100 Hz Clock) 
Polarimeter Left 

Polarimeter Right 

Not Used 

Not Used 

LBRF Coincidence 

LFRB Coincidence 

LFRF Coincidence 

LBRB Coincidence 

Polarimeter Left Accidentals 
Polarimeter Right Accidentals 

Not used 

Polarimeter Left+Right 
Polarimeter Left+Right (Prescaled) 
Polarimeter Leftt+Right Accidental 
Pulser (=Polarimeter Left+Right Prescaled) 
Not used 

LES CNaE *P' 1) 

LbaCNal sel) 

RB (Nal °P1) 

RE CNa TPL) 

Not used 

Nal LF 

Nal LB 

Nal RB 

Nal RF 

Not used 

Pulser (=Polarimeter Leftt+Right Prescaled) 
Event Fast 

Event Fast Coincidence Presented 
Pbasticakek 

Plastic) LE 

Plas tie RB 

Plastic LB 


* In cases where a parameter has two scaler numbers, the 
second number corresponds to an overflow scaler. 
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Word Number 


WOON UNF WNr OO 
I 
oO 


NPP RP RP PRP RPP Re 
OOMONADAUFWNHO 


Nw NY 
WN rr 


i) 
> 


25-26 
27-28 
29-30 
31=32 
33-34 
35-36 
37 
38-end 
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TABLE IIIf 


Event Record 


Parmeter 


-(Length of Event) 
Eline) (= ort a2.) 
LAM PATTERN 
SNe 

PRE Ho 

DGR (24> bits) 
TDC LO 

DG. eee: 

THC? 2 

LDCe Lo 

TDCe RO 

TDG RI 

TDC FR2 

TDC R3 

ADC LO 

ADC LL 

ADG~ L2 

ADC L3 

ADC L4 

ADC RO 

ADG RI 

ADC R2 

ADC R3 

ADC R4 
Scaler 
Scaler 
Scaler 
Scaler 
Scaler 

Sica ler 

MWPC 

MWPC 


Structure 


Comment 


Not used 
Not ‘used 


Lett* Start, Stop=-Nati ULF 
rf 
Nal LB 
R Plastic 

Ri ghia Start. Stop=NalskeF 
Ie AE 
Nal RB 
Beptas oie 

Plastic LF 

Nal LF 

Piactic LB 

Nal LB 

Spare-Not used 

Pease tse ck 

Nal RF 

Plastic, RB 

Nal RB 

Spare-Not used 

Ion Chamber 

Event Fast 

Master Gate 

Time 

Polarimeter Left 

Polarimeter Right 

Number? Of hits +°4 

Coordinates 


ipey . 307497184 sodnu¥ 
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were organized into runs. At the beginning (end) of each 
run, a header (trailer) block was written which included the 
contents of the scalers listed in Table Il. By subtracting 
the header block scalers from the corresponding trailer 
block scalers, the accumulated values for the run were 
obtained. For a given run, the angles of the detector booms 
were kept constant as was the polarization direction of the 
beam. Successive runs recorded on a single tape were 
separated by an end of file (EOF) mark, while the last run 
was followed by a double EOF mark. As well as processing 
the information and writing events on tape, the acquisition 
program recorded the scalers and kept a set of histograms 
for on-line analysis of the data. The spectra stored in 
memory are listed in Appendix II, along with their DCR 
gating requirements. 

During the experiment, the detectors were positioned at 
various angles by remotely controlling the booms. The 
pairs of detector angles for which coincidence data were 
acquired are shown in Table IV. Data for each row of angle 
pairs in this table were acquired simultaneously. CD> and 
*Oca data were acquired for the same angle pairs in order to 
allow further calibration of the Nal detectors using (p,2p) 
on /42C or deuterium. Data were taken with the CH, target to 
provide calibration points for the Nal detectors from free 
p-p scattering. Only the angle pairs 57°-30°, 30°-57° and 
44°-44°, from the various CH, angle pairs, provided 


calibration information. 
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TABLE IV 


Angle Pairs of Data Taken During Experiment 


*0Ca and CD, Targets 


LB-RF LF-RB LF-RF LB-RB 

54°-30° 30°-54° 30°-30° 54°-54° 
494—=35" 30°—54° 102-35. 49°-35° 
67°-25° 25°-67° 25" —235— 67°-67° 
42° -40° 30°-62° 30°-40° 52°-62° 
47°-29° 29°-47° 29°-29° 47°-47° 

CH, Target 

LB-RF LF-RB LF-RF LB-RB 

57°-30° 30°-575 30°-30° A da 
44°-25° 25°-44° 230-250 44°-44° 
63°-44° 44°-63° 44°-44° 63°-63° 


Using the +*0Ga tarvert, typically 3 to 5 tapes were 
acquired at each set of angle pairs and each polarization 
direction. Approximately two tapes (one third of a tape) of 
data were acquired for each set of angle pairs for the CD, 
(CH>) target. One tape contains roughly 150,000 events 


divided among the 4 detector coincidence combinations. 
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CHAPTER IV 


EVENT ANALYSIS 


4.1 Introduction 

The reduction of the information on the tapes to the 
final results can be divided into two phases. In the first 
phase, all events for a given beam polarization direction, 
detector coincidence combination, and angle pair were 
binned into two dimensional histograms containing prompt 
(realstrandoms) or random events. The histogram axes were 
T,tT>, (ESUM) and T,-T> (EDIF). T, and T, are the kinetic 
energies of the protons at the target. With appropriate 
corrections, these numbers were converted into cross 
sections. This chapter discusses the first phase of the 
analysis in considerable detail. Readers who desire only a 
general understanding of this phase are advised to read 
sections 4-1, 4.2, and 4.14. 

The quantity EDIF for three of the four detector 
coincidence combinations (LF-RB, LF-RF, LB-R3) was defined 
as Tieft-Tright- Since some of the angie pairs 
for the RF and LB detectors were the same as those for the 
LF and RB detectors, data for LB-RF was binned as Trr- 
TLB- With this convention and for similar angie 
pairs, a spin + run for LF-RB would be equivalent to a spin 
+ run for LB-RF - thus spin + (+) LB-RF data was treated as 
spin + (+) data. Due to different ESUM resolutions and 


accuracy of the energy calibrations between LB-RF and LF- 
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RB, events from these detector combinations were always 
binned in separate histograms. 

The data are organized in ESUM versus EDIF histograms 
because farpsinglef particles state dis techaracterizedy! bywea 
unique binding energy, BE. Leet, sist the tincident ‘enengy 
and T, is the energy of recoil for a (p,2p) reaction, ‘then 


As the recoil energy of the residual nucleus (T3) is quite 
small for +*°Ca(p,2p) Chessi¥than evOn55) MeVe foreerecoth 
momenta less ‘than 200 MeV/c), the cross» section for «a 
single*particlée stateigjshould* peak. tat atnearly fixed ‘value 
of ESUM independent of EDIF. Varying EDIF while keeping 
ESUM fixed changes the recoil momentum allowing the 
examination of the momentum dependence of a single particle 
state. As the’ cross section for a state changes ‘rather 
ealawhydasica fanction lofierecovl momentum; ithe sEDIF« bintsize 
51820) MeV 20 SEhewsESUMS bin fetize isichosen? to’ besl<0 MeVesince 
this is 1/3 to 1/4 of the summed energy resolution (FWHM) 
attained in this experiment. 

In the second phase of the analysis, which is 
discussed in the next chapter, cross sections for each of 
the ld3z,;5, 28,,9, and lds,» single particle states were 
@xteractea® By ietting Causst7antSpeakssotos ESUM ‘spectra 
for the different 20 MeV EDIF ranges. As there were data 
EGGStboth po barization t#directirons; the Yresulbtsipich )«the 


fitting were converted into unpolarized cross sections (59) 
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and analysing powers (Ay) for each single particle state 
as a function of EDIF and angle pair. 
4.2 Overview 

An outline of the first phase of the analysis follows. 
This portion of the data analysis was performed on a Data 
General Eclipse S200 computer. The tape reading software is 
described in (AN-78) while the remainder of the event 
analysis software was developed specifically for this 
experiment. For a given run, the data were subjected to a 
number of cuts and the surviving events were binned into 
eight ESUM versus EDIF histograms. By means of windows on 
the, left-right timingsi(start=left, stop=right plastic), 
events were binned either into a prompt event histogram 
(corresponding to the central peak of Figure III-6) or a 
random event histogram (corresponding to the left and right 
Deakess Otaeliourespiit-6 jie forjjeaeh),of Athermouriodetector 
coincidence combinations. The assumption was made that the 
energy distribution £ orn random coincidence events 
originating from two beam bursts was the same as that for 
random coincidence events arising from a single beam burst. 
Since. both. random. coincidence, peaks were included to 
increase statistical accuracy, the random event histogram 
contained twice the actual random background correction. 
The counts (divided by two) in the random event histogram 
were subtracted from the counts in the corresponding prompt 
event histogram before the cross sections were computed. 


GOrihtieria Ceuts)iwsed«tOr reject cevents «fromy theledata 
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were applied equally to the prompt and random data. Some 
of these cuts, although necessary, eliminate valid events 
and compensating corrections were necessary to avoid 
biasing the, final cross sections. For this reason, during 
the reduction of each run, statistics and correlations were 
acquired from which correction factors were calculated. 

A.a) Nal Detector Efficiencyeand Solid Angle Definition 

For a given particle energy, the efficiency of a Nal 
detector is nearly constant over its central region and 
drops off sharply near the outer circumference due to out- 
geattering 300 CA-77). Thegtradius of -the central plateau 
decreases ‘with increasing’ *particle energy. The energy 
dependence of the Nal efficiency in the plateau is given in 
able- V. For energies greater than 150 MeV, these 
efficiencies include the effect of inserting enough copper 
to degrade the energy to ~135 MeV. A correction factor 
Ponereting@or the product@o0G the Nal efficiencies for the 
detector energies of the two protons has been incorporated 
into the Galculation of theSeross sections. 

The MWPC coordinates were used to define circular 
windows centered on the faces of the four Nal detectors. 
The center of each of these windows was established by 
examining the distributions of WGatp.2p) “events ‘on x 
versus y histograms. The radius of these windows was 
chosen on the basis of histograms such as the one shown in 
Figure IV-l. This illustrates the radial distribution of 


(p,2p) events per unit area from the center of the detector 


48 


9 Oa) : ws » di ' 


* 4 , A y e : wz 
Wie »> ATR t ons a WwRE ¥ Pe It; , ; og vibe: pa det lene a ; 


lev. o ; $257 $7 3% & i200 Feo gelleen ‘gkeya 
Git. ee Puro i ety4 4 


‘ bizseadgaz y ioie To. geieootat od 
; ; ‘oe 


< ta601 £ segs P 7a 
i« i , y uj E t & .e } 44 
. = es Pee at if ; 2 i 70 
= bal : : 7 
a 
Py 7 = é - * es . * ' 7 = y _ , « , « ? = 
; on : . t6¥q. 2eeJea s ites i :e3569¢ 
, 7 
« _ 
i 4 - a = ‘ a - = 
ene , ‘i 56% is » Vlesato 270° eqeTh 
~ s - 
nee ‘ 1. ag rL TAS aie yD O22 
v _ 
’ 4 "> : 1 ee ite. €246 
WGaysos 3 to @4ne> 
. am 
2 + + 2 > si fi + 4 9% = :¥ 


iat i 0, 20 JewbASY Say Se pregetene 
; a ae | - 
bes — J ri ‘ snat« 7 i mel “a7 %% aut giend toSusse8 
: ear ‘. oa 
262 ero w elt In eeteiewlorios edt os 


“ 


wv: ibieb 4 bb oe (ee rieatiasoo hare re 


b3o 2 low eof en72 30) 1 ed eet eo seheahow evobels 


—— 
+ =—— 
a? ? 


uw betebi¢aree env Jewebite canis 1S fsao 20° vie AT 


2 ae Qveows (¢f,q)900"" Yo neki aT MANS wee tntana ad 
— 


afd pr we 


* Wobnty specu Sa wth oey kabel 


‘pita oo five ane biachiah 


, “ » fnabes' ‘os ; 


TABLE V 


Nal Efficiency versus Incident Proton Energy 


Energy (MeV) Efficiency 

(MeV) Eg 
10 0.998 

20 0.994 

30 0.987 

40 0.978 

50 0.968 

60 0.958 

70 0.947 

80 O93) 

89 0.931 
104 O93 16 
139 0.866 
147 0.844 
210 O.i27 
221 0.728 


Values for 10 to 80 MeV from (CE-69), 


for oo) to 221 MeV from (CA-/7). 


Uncertainty 


in EQ 


0.0050 
OS 0021 
0.0055 
0.0015 
0.0063 


0.0030 
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FIGURE IV=1 

Counts per unit area for each of the 4 Nal 
detectors as ‘a “tunctton) of Mdistancesmag 
wires, wire spacing = (2mm)* from) the™cenger 
of the detector-+ Counts are LF-“RB (30-5459) 
and . LB-RF (49°=35°) “events, with “a 40Ca 
target. Events falling outside a radius of 
26 wires (indicated by arrows) were rejected 
in the analysis. 
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- the abscissa has units of wires (where the wire spacing 
is 2 mm). For each of the detectors, the central region 
is flat (as expected) and falls off beyond 26 wires. Thus, 
a common window radius of 26 wires was chosen for each of 
the detectors. The resulting solid angles are LF=4./2msr, 
RF=4.73msr, LB=6.36msr, and RB=6.43msr. Events outside 
these windows were rejected. The front (back) detectors 
subtended ~12.6° (~4£2.2°) at “the target. 
4.4 Multiwire Efficiency 

Wires away from the edges of the multiwire planes are 
quite efficient (typically >99.5%). However, any event for 
which one (or more) of the four wire planes fails to 
register a hit must be rejected. For a certain fraction of 
these events, this missing coordinate would be within the 
window. Consequently, a correction consisting of the 
product ‘of © the’ Seftficiencies “of "ithe. *four “participating 
multiwire planes is included in the calculation of the 
cross sections. For the x and y planes on one arm, the 


efficiencies are given by 


XMISS YMISS 
GEV£2/) MWEFF = 1 YaIT . peeps es = 1 tiir 
where XHIT,YHIT are the number of events within 26 wires 


of the center in X (Y), 
XMISS,YMISS are the numbers of missings in X (Y) 
whose corresponding Y (X) coordinates were within 
26 wires of the center. 


Typically, the product of the efficiencies of the four 
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planes was ~0.99. 
4.5 Pile-Ups and Multiple Hits 

In the event analysis, events with a pile-up DCR bit 
set for either detector arm were rejected. These 
represented ~13% of the total number of events. 
Approximately 0.5% of the events had. pile-ups in both arms. 
In adjusting the cross sections for the pile-up cut, each 
pile-up was counted as representing a single (p,2p) event. 
Any double (p,2p) events rejected by the pile-up criteria 
were automatically accounted for in the computer dead-time 
correction. 

When a proton passes through a MWPC, it is most likely 
that a single wire will be struck, although it is possible 
that two or more wires (not necessarily consecutive) will 
be struck. It is also possible that two particles may pass 
Ennouch Stes iplane-awlthin the -resolving time,  tyy, of 
the plane (~100ns) causing two distinct sets of struck 
wires. Furthermore, two particles passing through the 
telescope sufficiently close in time will fail to be 
resolved by the pile-up monitor and thus the pile-up DCR 
bit will not be set. 

The span of a struck plane is defined as the number of 
wires from the leftmost (or topmost) struck wire to the 
rightmost (or bottommost) struck wire. Thus, a single 
struck wire has a span of one, two struck wires separated 
by a gap of two wires has a span of four and so on. If the 


assumption is made that only a single particle passes 
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through each wire plane per event, then a reasonable 
definition of the position of this particle is the sum of 
struck wire coordinates divided by the number of struck 
wires (centroid). However, for large spans, the concern 
arises whether two particles have in fact passed through 
the plane or whether a single particle has produced a large 
span (for example, by ejecting a é-ray to one side). In 
either case, as the span increases, confidence that the 
calculated centroid represents the actual centroid of the 
particle(s) rapidly erodes as there is no reason to assume 
that the span is symmetrical about the position of the 
particle(s). 

The distribution of events for the y plane of the LF 
MWPC is plotted as a function of span in Figure IV-2. This 
distribution is typical of all the wire planes and it was 
decided to reject events for which any one of the 4 planes 
recorded a span greater than 4 wires. Such events are 
referred to as multiple hits. After pile-up events were 
rejected, 1224 of the remaining events were found to be 
multiple hits including 0.5% which had spans greater than 
four on all four wire planes (quadra-multiple events). No 
account was taken for the fraction of the quadra-multiple 
events which may have been double (p,2p) events and all 
multiple hit events were treated as single (p,2p) events 


with the final cross sections being corrected accordingly. 
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MECURE MLV = 2 
Distribution of events (expressed as_ the 
percentage of total) passing through y plane 
Cif atee2 ama Seca eeeuuc co of span. Events 
with spans greater than 4 (indicated by the 
arrow) were considered multiple hits. 
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4.6 Elimination of Free p-p Beare 

The calcium target used in the experiment contained 
substantial amounts of hydrogen (H/*°a ~ 0.23 in atoms) as 
discussed in section 4.13. All the angle pairs for LB-RF 
and LF-RB (with the exception of 47°-29° and 29°-47°) 
overlapped with free p-p scattering angles. The relatively 
large free p-p cross section (~10.6 mb/sr at 200 MeV at 
43.6°) resulted in the presence of considerable numbers of 
free p-p events in the data. The peak for free p-p events 
is well separated from the (p,2p) data as shown in Figure 
IV-3 where Tyeft versus Trignt is plotted for 
30°-54° (LF-RB). However, energy losses’~ suffered ro 
elastic protons undergoing nuclear reactions in the Nal 
erystals result in two p-p reaction tails which extend into 
the (p,2p) region as can be seen in the figure. Free p-p 
events were eliminated by means of coplanarity and opening 
angle cuts applied to the data based on the distribution of 


Ax versus Ay where 


Gives.) 2) Ox. = etm entice ; 
Cy BP ban ) 
ete lett, “right ° 


2 


NON ely ee COLrTea pond std, asthe ‘opening ;angles. in the 
horizontal and vertical planes respectively. Such. 7a 
distribution for events from the free p-p peak (isolated 
from the *Ca(p,2p) events by means of applying suitable 


windows to the Nal signals from both arms) for 30°-54° is 
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PoUGU RE eee Via 
Scatter plot (logarithmic density scale): for 
Ti, versus Tpp- Counts are LF-RB (30°-54°) 


events with a +a target. Tiel Cp, 2 p)peloes 
as well as the free p-p peak and reaction 
tails are labelled. The reaction tail on 
the low) enerey <side LGB ets son lyemba pew 
visible. 
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shown in Figure IV-4. Eliminating all events which fell 
within a rectangular window surrounding this distribution 
(shown in Figure IV-4) resulted in the Tjeft versus 
Tright distribution shown in Figure IV-5 (which is to 
be compared with Figure IV-3). This method, which we refer 
to as the elastic cut, required a different window in Ax 
versus Ay for each angle pair of LB-RF and LF-RB. St 
eliminated >96% of the free p-p events. 

However, it also excluded a fraction of the (p,2p) 
events. To estimate what fraction of (p,2p) events were 
being rejected, the elastic cuts were applied to runs not 
containing free p-p events (47°-29° and 29°-47°). By 
counting the numbers of events eliminated, we obtained a 
correction factor for each elastic cut. This correction 
factor was incorporated in the calculation of the cross 
sections. The size of the correction varied with the Ax 
versus Ay window size and ranged from 1% to 19% due to 
varying degrees of Rafe with free p-p angles. 

4.7 Nal Calibration and Energy Resolution Optimization 

it Sorder “to “bin ~events .in the... ESUM versus. EDIF 
histograms, pulse heights from the left and right Nal 
detectors were converted into energies. Lo.sdo. sthiss:;-) 4 
mumber of factors had to be taken into account. The 
relationship between Nal pulse height and energy was 
nonlinear. In addition, the Nal ADC pedestals, zero 
offsets due to a de level added to each Nal signal, were 


not negligible. Since the accuracy in determining the 
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FIGURE IV-4 

Scatter plot (logarithmic density scalejmroc 
Ax versus y/2%) (proportional to ~openiag 
angles in the horizontal and vertical planes 
respectively). Counts» are 9 GhokBe30.—o4ee) 
event easwlth amt Ca target. The free p-p 
events (selected by software windows) and 
the few remainine ssp, 2p) ee Vien: saumeoimae 
labelled. Events Palbine® “within the 
rectangular ‘region ~were ‘rejected sein e=ethe 
analysis. 
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BEGURES IVa 

Scatter plot (logarithmic density scale) mon 
Trp versus Ton. CONnditions=are tne ss amemas 
for FiguressiV—3) bute witheablmevemesmrat aa 
within the rectangular window shown in 
Figure IV-4 rejected. Although >96%4 of the 
free p-p events have been eliminated, the 
p-p peak is still visible because mote cve 
logarithmic scale. 
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cross sections of the single particle states depended 
strongly upon the degree to which the peaks were resolved, 
it was worthwhile to invest considerable effort in removing 
nonlinearities and minimizing effects that degraded the 
energy resolution. 

Thescfrees p=p scattering in ithe “LB=RFY and 4LF-RB +*0¢a 
runs and in the CH, runs provided a number of well-defined 
energy calibration points. Free p-p events were isolated 
from the *0Ga(p,2p) data by means of Nal pulse height 
windows and, in addition, the left-hand-side protons were 
required to pass through a 6 wire wide window in the x 
direction. Since scattering angle and energy are 
correlated for free p-p scattering, this ~0.6° window 
served to define precisely the energies of the protons at 
the target. The position of each x-plane window was chosen 
and a restriction was applied to the y-coordinates in order 
to guarantee that the protons in either arm passed within 
the circular region of constant Nal efficiency. This 
minimized the reaction tails, making the precise 
determination of the peak positions easier. 

Protons emerging from the target passed through a 
considerable amount of material (the target itself, the 
scattering chamber window, air, the MWPC's, the plastic 
scintillators, and caps on the fronts of the Nal detectors) 
losing some of their energy before being stopped in the Nal 
detectors. For calibration purposes, given the energies of 


the free p-p protons after the interaction in the target, 
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their final values at the Nal detectors were required. 
Once the detectors were calibrated, the energy deposited in 
the Nal detector was used to obtain the scattered energy at 
the target. For both purposes, a table of values 
representing the energy-degrading materials between the 
target and the Nal detectors was necessary. A suitable 
table was compiled and is shown as Table VI. The energy 
losses experienced by protons passing through these 
materials were obtained from tables of stopping powers (SE- 
67). Energy losses suffered by 146 (54) MeV protons at the 
target reaching the front (back) Nal detectors were ~8.1l 
(~13.4) Mev. 

For each CH, run and for the p-p elastic events in 
each +*%a run, the average value of the NaI pulse height 
for the free p-p peak was obtained. In addition, the 
average LED and pedestal pulse heights for each of the Nal 
detectors for these runs were obtained. A plot of average 
pedestal pulse height versus run number for each detector 
is shown in Figure IV-6. The common fluctuations among the 
four detectors are believed to be due to a combination of 
temperature effects and electronic noise while the origin 
of the similarity between the left or right pedestals is 
believed to originate in the ADC's. The fluctuations over 
the course of the experiment are on the order of 5 channels 
(#1.10 to 1.42 MeV depending on the detector) while the 
resolution (FWHM) for a given run ranged from 4 ‘to 7 


channels among the various detectors. 
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TABLE VI 


Materials* for dE/dx Energy Loss Calculations 


Region Thickness *Density Material Comment 
(gm/cm*) 
LForRF LBorRB 
1 Dependent on "lGa, Target Energy Losses. 
angle and target CD5,CH>. 
2 0.134 OE V2 Air Scattering Chamber 


window to front wire 
plane window. 


3 0.033 0.033 Copper MWPC wires and gas. 
4 0.010 0.015 Air Back wire plane 
window and air gap. 
5 0.835 0.598 CH Plastie Scintillator. 
6** 5.6.00 ab F = Copper Degrader. 
5.644 RF = 
7 0.583 0.250 CH» Material in front 


of Nal crystal. 


* Note: Some materials in table are the stopping power 
equivalents of the materials they represent. 


**Ueed only for LF, RF = 25*.. 


TABLE VII 


Materials for dE/dx Energy Loss Calculations 
Adjusted for Nonlinearity Correction 


Region Thickness ‘Density (gm/cm%) Material 
1-4 Same as in Table VIII 
LF RF LB RB 
5 0.500(0.835*) 0.835 0.460 0.295 CH 
6% SL 0 5-644 = = Copper 
7 - 0.183 - - Cy 


*4Used only. for ‘LF, REY= 25°. 
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FIGURE IV-6 
Nal pedestal pulse height centroids for each 


detector as a function of run number. Data 
shown .are for "Ga runs. white gaps 
correspond 9to) acquisition ot 96GH (O tee UD, 


data. 
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LED pulse height (corrected for the pedestal) versus 
run number is plotted in Figure IV-/7. The large shifts in 
gain (~20 to 45 channels), which are correlated among the 4 
detectors, are very probably due to stability problems in 
the common pulser input to the LED's. The resolution of 
the LED peaks (FWHM) for a given run ranged from 11 to 16 
channels among the four detectors. Finally, the average 
pacisenonedgeit -of cthe sfiree “pp events ‘Ccormrnected for 
pedestal) versus run number is plotted in Figure IV-8. The 
ESUM resolution (FWHM) varied from 20 to 24 channels. Runs 
with the same detector angles are grouped together in the 
Puscure “since ~the average pulse height for free p=p 
scattering varies with angle (data for Figures IV-6 and IV- 
7 were organized similarly to facilitate comparison). For 
data corresponding to a given set of detector angles, the 
variations are probably due to a combination of ADC 
fluctuations and individual detector base noise. The 
changes in the LED pulse heights are substantially 
larger than, and not obviously correlated with, those of 
Bue free p-p events. EOr- this “reason, “the -LED- pulse 
heights were not used in the calibration of the Nal 
detectors. If the signal as output directly from the 
pulser driving the LED's had been recorded on an event-by- 
event basis, it might have been possible to remove these 
variations and the LED pulse heights would then have been 
quite Useful in the energy callbration. 


The method of calibration was as follows. initially, 
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FIGURE IV-/7 
Nal LED= pulse  heleht) (minus pedectaa® 
centroids for, each detector asma tinction mon 
run number. Condittons actor igure Veo. 
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FIGURE IV-8 

Free p-p Nal pulse height (minus pedestal) 
centroids for each detector as a function of 
run number. Data are from scattering from 
hydrogen contamination in + Ga target. For 
each set of data, the detector angliewand= the 
difference between the measured pulse 
heights and :the plotted values are given 
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pvewmenerteyociac ) the "Nel sdetector  (l.41), ~<as determined 
from the pulse height, was assumed to be related to the 
pulse height of the proton (PHp) and pedestal 


(PHped) in a linear fashion, 
Give T= x(iycPHoCi)-eH 4 G4)] 


where g(i) is the gain for the ith detector. 


The gains were obtained from the average pulse height for 
free p-p events, pret (iy, and pedestals, PH ag (ht) from the 
P 


44°-44° (LF-RF, LB-RB) CH, runs: 


ee) 
OD) BG) te ee 
[PH (4)-PH (4) 1 


where Teaintt? is the energy at the ith Nal detector as 
determined from kinematics at the.) target 
x (corresponding to the angle of scatter) and 
corrected for energy losses. 

With this set of gains, a value of energy based upon 
the average free p-p pulse height and pedestal for each of 
the detectors for each angle pair of *0Ga data containing 
free p-p scattering was computed using equation IV.4. The 
result of these calculations is presented in a plot of 


versus shown in Figure IV-9. The 


es eres 5 Pkin 
departures from zero at high and low energies are a result 
of the nonlinearity of the system. A nonlinearity 
correction, of the necessary energy dependence, was 


introduced by adjusting the amount of material in front of 
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FLGURE LV 9 
Kinematically determined energy at Nal 
minus the energy determined from the Nal 
free p-p pulse heights Gee Tat? £Ox 
each of the 4 detectors as a function of 
T,..+- Pulse heights are from CH. as well as 
kin 2 
*“[Ga runs. 


HPEGURES LV LO 
Same as in,;*Figure 1V-9 Mafter Snonlinearii, 
correction has been applied. 


Tain ~ Toac (MeV ) 


OLF *LB 
<RF *RB 


0 20 40 60 80 wee ae eee 


¢ 7 y 
ah ¥ 
ee CY 


ee 


4 


| ” 
» P 
vo ' 
Pe oe he 4 
oA » Ww 
= * = / 
# _ s 
i 
. ! 
j y= 
€. » r 
1 
od 4 
+ a 
e*: 
) 
-_ = 
« 
>) ’ 
@ 
, 
¥ 
If 
; 


the Nal detectors used in the energy loss calculations. 
The adjusted values are given in Table VII. 


With the same data, the gains were recalculated using 


ad j ad j 


: owt 
id EeRet Sethe 


equation IV.5 where T was replaced by T 


kin 
energy at the Nal obtained from free p-p kinematics at the 
target and corrected for energy losses using the adjusted 
densities given in Table VII. The process used to obtain 
the points for Figure IV-9 was repeated and the new energy 
calibration gives the results shown in Figure IV-10. The 
non-linearities are now seen to be largely removed. 
Thereafter, energies at the target were obtained from Nal 
pulse heights and pedestals by ‘means of equation IV.4 
(using the adjusted gains) and by performing energy loss 
calculations with the values from Table VII. 

With this “procedure, the CD, target data were binned 
in ESUM versus EDIF histograms. The resulting positions in 
ESUM of the !2c(p,2p) (BE=15.95 MeV) and *H(p,2p) (BE=2.225 
MeV) loci typically ranged from 3 MeV above to 5 MeV below 
the expected position of the loci (based on kinematics and 
taking into account the recoil energies of the 11g and 
neutron). The expected positions in ESUM of the loci for 
12¢(p,2p) and 2H(p,2p) ranged from ~182 to ~184 MeV and 
from ~190 to ~19/7 MeV respectively. For a given energy 
difference and angle pair, the difference between the 
position of the data and the expected position was the same 
for the 12¢(p,2p) and 2H(p,2p) data. This, along with the 


fact that the size of the bins in EDIF (20 MeV) was large 
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compared with the error in ESUM, led us to regard any 
furtner eadijustment ijof the energy calibration as 
unnecessary. The remaining discrepancy was regarded as a 
linear shift in ESUM which was treated in the fitting of 
the data. 

A further refinement of the pulse height to energy 
conversion process involved the pedestal values. The 
pedestal for each detector was continually updated by the 
average of the last 100 pedestal values for that detector. 
This resulted in an improvement of 0.2 to 0.3 MeV in the 
ESUM resolution. The ESUM resolution varied with the angle 
pair and detector combination and ranged from 3.2 to 4.4 
MeV. The variation with detector combination reflects 
individual detector differences. The variation with angle 
pair is possibly caused by time-varying temperature effects 
on the detector electronics. 

4.8 Reaction Tail Subtraction 

Sate es eee the Nal(T&) crystals produce tails from 
every point along a (p,2p) locus similar to those shown for 
Meee ap pesscatteringwin Figure .[V—11-. To ensure that the 
ESUM resolution of the individual valence states was not 
being seriously degraded and that the heights of their 
peagre, G10 pnOCe nave -alsienificant background dué-,to’.the 
reaction tails of higher energy protons, an examination of 
the effect of reaction tail subtraction was performed. 
The shape of the reaction tails was taken from p-p spectra 


from the CH, runs and an attempt to strip away the- tails 
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ELGURE 1V—11 
Scatter plot (logarithmic density scale)s tor 


Trg versus Tp,- Counts are LB-RB (44°-44°) 
events with a CH, target. The free Poe peak 
and reaction tails as well as the !?! COP 2 pp 
locus are labelled. 
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from the (p,2p) data was made. For a given angle pair and 
spin direction, (p,2p) events were binned in a ESUM versus 
EDIF histogram with bin sizes of 1 MeV by 1 MeV. For each 
bin of the histogram, the reaction tail contributions due 
to all higher energy protons were subtracted before 
proceeding to the lower energy bins. Once the subtraction 
process was completed, the remaining counts were rebinned 
in an ESUM versus EDIF histogram with bin sizes of 1 MeV by 
20 MeV and compared with the uncorrected histogran. The 
reaction tail subtraction made no measurable improvement in 
the ESUM resolution, and an upper limit for the decrease in 
the counts for the ld3,5, 28),9, and ldco,» region was <22. 
Thus, the effect of tail subtraction was small and was not 
included in the final results. 
4.9 Removal of Deuterons 

Events from the (p,pd) reaction were acquired 
concurrently with the (p,2p) events. To eliminate them, a 
Geatiit yoriproporcaonal to. “the, masses -of “the detected 
particles was calculated. The range, R, of a particle of 
energy T and mass M in a given material can be approximated 
by (SE-64) 


(Iv.6) R(Z) = K-(5)" , 


where K,n are constants, and n*1l.8. 
From an expression relating the difference in range (AR) 


for two particles differing in energy by an amount AT, 
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and using equation IV.6, the following expression for the 
mass may be derived: 


(Iv.8) Me [crt ar)t-¢r)" ]t/® : 


Here AT and T represent the energy deposited in the plastic 
detector and the emerging particle energy respectively. 
Thus, plastic scintillator pulse height had to be converted 
to energy for each of the detectors. Pedestals and free p- 
p scattering pulse heights for the plastic scintillators 
from CH» runs were used to establish a linear relationship 
between pulse height and energy. The energy of the 
particle venerging from ‘the plastic -scintillator was 
obtained from the energy deposited in the Nal detector plus 
the energy lost in travelling between the plastic and the 
stopping counter. 

In order to resolve protons and deuterons on the 
basis of mass, it was necessary to correct the pulse 
heights in the plastic scintillators for an observed 
dependence on the vertical position of the particle in the 
detector. This dependence was studied by examining free p- 
p events from the CH, runs and the results are plotted in 
Figure IV-12 as a function of vertical distance (in wires) 
from the center for each of the plastics. 


This behavior probably arises from a number of 
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EPGURE LV s1eZ 

Relative. plastic pulse heights for each 
detector as a function of vertical distance 
from ‘ithe ‘center “of the “event-detinime 
window. Data are from free p-p scattering 
from the CH») target. The postition sof athe 
light guide and the center of the window are 
indicated. 
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effects. The Siight sundersoes ‘attenuation: in ‘the 
scintillator, and the solid angle acceptance of the light 
guide decreases with increasing distance from the guide. 
In addition, light signals undergo more reflections from 
the boundary of the scintillator as they propagate toward 
the light guide the higher they: are generated in the 
detector, and with each reflection there is some loss of 
signal strength due to partial internal reflection. This 
iucerpretation= correctly predicts ‘that (the “thick: front 
plastics should experience a smaller effect than the thin 
back plastics (which should require more reflections). 
After correction for this position dependence, the 
protons were sufficiently well resolved to allow 
elimination of the deuterons. Figure IV-13 shows’ the 
distribution of events as a function of the mass parameter 
for the RB detector for 30°-54° ‘*%Ca data. The deuteron 
peak is ‘to the right of the larger proton peak. Mass 
parameter windows for protons were established for each 
detector in a (p,2p) coincidence and events with either 
particle mass outside these windows were rejected. This 
eliminated varying fractions of events from ~9% (25°-67°) 
to. “262 (67. -67~'). Generally, the angle pairs for which 
typical (p,pd) recoil momenta were smaller than typical 
(p,2p) recoil momenta had the greatest number of events 


containing deuterons. 
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REGURESSLV— ts 
Counts as a function of the mass parameter 


for the RB detector. Counts are LF-RB (30°- 
54°) events witheaws Ga target. Proton and 


deuteron peaks are labelled. 
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4.10 Decelerated Beam Events 

The fraction of decelerated protons in the incident 
beam was large (~12%-78%) in 7 of the 39 *%Ca runs. In the 
remaining 32 runs, the decelerated beam component was <2. 
As mentioned previously, decelerated beam events could be 
resolved from normal beam events by means of the timing 
versus rf TDC information. Using the same procedure as in 
the Nal calibration to isolate free p-p events from *%Ga 
data, the energy spectra for decelerated and accelerated 
beam events were studied separately. It was concluded that 
the incident energy of the decelerated component of the 
beam and the ESUM resolution were indistinguishable from 
those for the normal beam. Also, the beam polarizations 
were similar to those of other runs. For these reasons, 
decelerated beam events were not separated or treated 
differently from normal events in the subsequent analysis. 
4.11 Dead Time Correction 

The loss of events due to computer and electronic 
dead-time must be compensated for in the calculation of the 
cross sections. Two (correlated) measures of the dead- 
time, the first (a) including only computer dead-time, the 
second (b) inchudings that a:fiiriom rtboth sources, were 
calculated from 

Events presented to MASTER GATE 


al s = ee tie. S| Ue ee 
(IV.9) DTC. Events written on tape ; 


prc. = Pulser events presented to system 
b Pulser events written on tape i 
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These measures were usually quite consistent; (a) was 
typically 20-62 to 1782" iarger than ‘(Cb)’: Prescaling of 
polarimeter counts caused the interval between subsequent 
pulser triggers® to be very regular. Thus, pulser events 
could cause dead time for non-pulser data but not 
subsequent pulser events. Ths), along weth= che*sfact*® that 
the electronic dead-time is expected to be negligible, 
explains the sense of the difference between the two dead 
times. 

The value determined by method (a) was’ used. In 
instances when there was serious disagreement between the 
two methods, the problem was believed to be due to a 
hardware malfunction in the pulser electronics. In such 
cases, the ratio of non-pulser events written on tape to 
non-pulser events presented was taken as the dead time 
eernrectiou. "Typically, this correction ranged from’ 1:05*to 
dit tor behe *°catirunas: 

4.12 Beam Normalization 

The number of incident protons during a run could be 
measured either by means of the polarimeter counts (NI,) 
OortAaon chamber*counts (NI},) “using *equations' IT'2**or 11.3 
respectively. The free p-p scattering cross sections at 
~43.2° from LF-RF and LB-RB CH, data with a polarization of 
~0.73 (+t) were computed. Using NIg, the values obtained 
were consistent to within 0.2% and were 5.0% lower than the 
value expected from phase shift calculations. When NIy 


was used, for LF-RF and LB-RB the results were 324 and 20% 
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lower than the expected values. 

Throughout the experiment, the ratio, of NIg to NIp 
varied from 0.85 to 0.95 while this ratio dropped to values 
as LoWnas 20.5 eto, 0. 7/9 shores, .eamalleinumberyoft suns. »ringthe 
latter cases, the event rates presented to the acquisition 
system were much more stable as a function of NI, than as 
a function of NIp. The variation in the ion chamber gain 
was believed to be due to beam movement, changes in beam 
halo, small changes in the gas mixture, the window bulging 
under increased gas pressure or variable atmospheric 
pressure, and, probably most importantly, impurities (air) 
leaking into the oe Because of its superior stability, 
the polarimeter was used for beam normalization. 

4.13 Target Contaminants 

Pure calcium is extremely reactive with many elements 
and compounds, including oxygen, making it notoriously 
difficult to protect from contaminants. Between 
experiments, the 20 Ca target was kept immersed in a mineral 
oil bath (primarily consisting of CH,) to protect it from 
aire After it was cleaned with acetone (CH ,COCH3), the 
target was mounted in the target ladder and lowered into 
the scattering chamber which was then pumped to a rough 
vacuum (7°90 millitorr = 12 Pascals). 

The amount, of. hydrogen inthe ‘target, at. the »time of 
the experiment, expressed as the ratio of hydrogen to 
calcium atoms, was O23 3002. This was determined by 


isolating the free p-p events from 40ca runs in the manner 
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previously outlined in section 4.7. Then, from known 
values for the free p-p cross section (from phase shifts), 
the amount of hydrogen could be deduced. 

The source of hydrogen in the target was probably 
Ca(OH)> and residual mineral oil not removed in the 
cleaning. As previlously mentioned, the free p-p peak is 
well separated from the *0Ca(p,2p) data and the reaction 
tails were almost entirely removed by means of the elastic 
cuts. However, while the free p-p events posed no problem 
in the analysis of the *0Ga Cp, 2p.) data, other contaminants 
such as oxygen (Ca(0H)>, CaO), nitrogen (CaN), and carbon 
(CaCO3) have (p,2p) loci which overlap those from 40Ga to 
varying degrees. Previously measured binding energies for 
single particle valence states for Can 1€, and@ °C are 
listed in Table VIII. Because of the 3.2 to 4.4 MeV summed 
energy resolution and the resulting overlap of single 
particle states from different nuclei, it is impossible to 
distinguish between the + 0Ca(p,2p) peaks and those of the 
contaminants. Consequently, other analysis methods were 
employed sing an attempt to ~ascertain the amount of 
contamination. Since these analysis methods were performed 
after various periods of time had elapsed since _ the 
experiment (performed in June 1977), the resulting values 
serve only as upper limits. Due to a number of 
difficulties, the methods were not equally sensitive to the 
various contaminants and the following represents the best 


estimate for each of the contaminants. 
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TABLE Vilt 


Proton Binding Energies for SoCas 160, 12¢ 
and Spectroscopic Factors for *0Ca Valence States 


Nucleus State Binding Excitation Sue er oecop = 
Energy Energy Factors** 
(MeV) (Mev)* 
e'ca Hees 8.33 0 ok 
28/2 10.85 Died 2 L165 
lds yo 13.60 Cer ay f 0.91 
13.94 5.61 0.64 
14.67 6.34 BTS Ae 
14.83 6.50 O10 
eye 6.77 0.10 
L529 6.96 0.138 
Bie ieee v) 0.10 
15.76 7-43 0.47 
1Ote ed (ees 0.10 
16.50 Sia / 0.35 
16.76 8.43 0.24 
16.88 pre Be Oras 
LLY Sees | 8.90 0.10 
Eie43 BeEO 0.10 
18.08 9.75 0.10 
4.96 
a0 ae Pe Pea127 3.797 
lp3 2 18.44 10.01 
==C Be ice Coit S 7.62 


*jaRelative, togld5;,5 state,.of Ca. 
** Spectroscopic factors and binding energies for *0Ca from 
(DO-76). 
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feeevalwe: e100 tne anterogen “to .calcium iératio of 
0.091+.004 was obtained (in August 1979) from a 200 MeV 
(p,d) experiment involving the +*0ca target and a 90.2 


2 kapton target (C 55H) —N 205, density=1.42gm/cm °) 


ng/cm 
which served as a nitrogen target. The scattered deuterons 
were detected by the TRIUMF Medium Resolution Proton 
Spectrometer (MRS) at 22.0°. TA the experiment, the peak 
for l4n(p,d) in the deuteron energy spectrum was identified 
by means of the Q-value for the reaction (*-8.3 MeV). The 
Kapton target was used to determine the cross section for 
this process and this in turn was used to determine the 
amount of !4*N in the +a target. Unfortunately, the 
peaks corresponding to Bien a) and 169(p,d) were not 
resolved from e’Capd)) peaks and no reasonably accurate 
value for carbon or oxygen could be ex npecteds 

A value for the carbon to calcium ratio of 0.054+0.009 
Was ms sodDtaumed Cin Aupust ~ 1979.) by means of gas 
chromatography. The method involved grinding a small piece 
of the target (~10 mg) into a powder which was then fed to 
a gas chromatograph. The sample was heated to a 
temperature of ~1100°C and the expelled gases were carried 
along by a jet of helium to molecular scrubbers which 
trapped = 2,0, COs, He, No, or. 05. Even though CaCO, 
dissociates at ~900°C, CaO would probably remain as a solid 
(dissociation point >2850°C). Also the dissociation point 


Cree caN, wets? 1192 Clocaada thus this. techni que was 


inappropriate... for determining _the oxygen ofr . nitrogen 
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contamination. 

The amount of oxygen was measured by scattering 18 MeV 
protons from the 40Ca target (University of Washington, 
Seattle, September 1979). Elastically and inelastically 
scattered protons were stopped at 110° in a solid state 
silicon detector. The ratio of !° to *%%a (RAT) was 


obtained from 


40 
G 
No Fao 
an 
GCLV..-k:0) RAT = ° : 
169 
N do 
G sei 
‘: af. 
where N N were the number of counts in the 40ca and 


Gan 0 
169 elastic peaks, and the elastic scattering cross 
sections, — for *%a and and 1°0 were obtained from (SH- 
f45VA-71).- The elastic peaks for 12¢ and ‘!4N were not 
resolved in this experiment. 

The * a target was briefly exposed to air in the 
process of mounting it in the scattering chamber despite 
attempts to isolate it by means of keeping it in a dry 
argon atmosphere. Thus, the value obtained for the ratio 
of 19 to +®ca (0.52640.055) represents a very pessimistic 
upper limit. 

In summary, the level of contamination for 14N and 12¢ 
appears low enough to be unimportant. However, the 
situation for !°®0 is less clear on the basis of these 


measurements. Fortunately, an indirect method of measuring 


the oxygen was found by means of incorporating previously 
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measured 169(p,2p) eross sections into the fitting of the 

*“Yoa(p,2p) data. Titseuemethod’ is tulby outlined in 

Chapter V and yielded a value of 0.324140.078. Furthermore, 

asp vdetailed in that ‘chapter, the final "Van. 2p) results 

were fully corrected for the presence of ay 

4.14 Event Analysis Summary and Cross Section Calculations 

The cuts and windows applied to the events from each 

*Oca run in the order in which they were performed are 

summarized below. 

i) Pulser events were rejected. 

ii) The detector coincidence combination was determined. 
Subsequent analysis focussed on information from the 
two detector arms in coincidence. 

BBY Events not passing through the circular multiwire 
windows were rejected. 

iv) The elastic cuts were applied to those runs and 
angle pairs ‘which-«contained free p-p events. This 
eliminated most free p-p events. 

v) the, -lett—richt “timing “(start=left,  stop=right 
plastic) was examined and events were flagged either 
as a random (peak 1 ‘or’ 3 of Figure I11I-6) or a prompt 
(peak 2). All subsequent cuts were applied in an 
identical fashion to both types of events. 

we) Events tor which pile-up occurred in elther detector 
arm were rejected. 

vii) Multiple hit events (events with spans >4 wires) were 


rejected. 
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viii) The energies deposited in the plastic and Nal 
detectors were computed. From these the mass 
parameter for each particle was calculated. If 
either particle fell outside the mass parameter 
window for protons, the event was rejected. 

ix) The appropriate random or prompt energy histogram bin 
was incremented. 

The bins of the energy histogram were: 100.5, 01254 
al:ews..5 209.5 for ESUM (where bin 100.5 covers the region 

1005101 ewe tetrcl <j; —160., p—b6 02,7 e2.2.0, 200 for EDIF (where 

bin -180 covers the region» [-190,-170) MeV: etc.). After 

thes aniarmation S@feom wmitliveiruns having. -a given ‘spin 
direction, nie Heentat: combination, and angle pair were 
combined, the cross sections were calculated. As the 

incident beam was only ~/0%Z polarized, the resulting spin t 

and ¥ cross sections are referred to as partially polarized 

cross sections. 
For a bin with energy coordinates ESUM=ES, EDIF=ED, 
bie epattcial hy potariwedser cross” "sectionga ands ennor . were 


computed as 


dito NE-NRE 
uaa) dESUM*dEDIF*d®,*dQ, — ( 2 anes 
Ad to = (NE+NRE ji/2 is 
dESUM *dEDIF *d 8) *d 2%, 4 ; 


where f is given by 


DTC *PUC *MHC *ELC 


(Iv.12) £& OESUM*AEDIF -dQ)°d25°N, °N| *MWEFF *NalEFF : 
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where 


where 


NE is the number of events in the prompt histogram 
bin, 

NRE is the number of events in the corresponding 
randoms bin (includes 2 randoms peaks), 

DTC is the computer dead time correction, 

PUC is the pile-up correction, 

MHC is multiple hit correction, 

ELC is the correction factor accounting for the 
elimination of (p,2p) events in the course of 
eliminating free p-p events; it applies only to 
those detector combinations and angle pairs with 
free p-p events present, 

AESUM (AEDIF) is the dimension of the bins in ESUM 
CEDILED. 7.1) 2.0;) Mev, 

d2,,d%, are the solid angles subtended by the MWPC 
windows, 

Nie Leet hewnumbers sol protons. sincident, on ‘the 
target (based on polarimeter counts), 

N; is ‘the number of *%a atoms per unit (area 
fugthe target «| 7644 20.97) x10 + cna, 

MWEFF is the product of the efficiencies for the 
four multiwire planes for the detector coincidence 
being considered, 

NalEFF = NalEFF(T',)*NalEFF(T',) , 

NalEFF(T) is the Nal efficiency at a proton energy 
| (from Table V), 


T',,T'5 are the energies of the protons at the Nal 
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detectors corresponding to energies T,, T, at the 
target, 

ES=Tj7T 5, 

ED=T ,}-To5- 

The expression for the, error in the cross section 
contains only contributions from counting statistics and is 
derived according to equation AIII.1 in Appendix III. A 
summary of the errors in the factor f and its constituents 


is wciven in Table 1X. 
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TABLE IX 


Errors in f and its Constituents* 


Typical Systematic or 


Statistical : Non-Statistical 
Uncertainties** 


Uncertainties** 


Af ea deer ere 2 4 LiDe See 
Component of f 

DTC A OM a4 0.6% 
PUC On. 6tee 40 38-4 

MHC OSbiae Gein. 92 

ELC 0.61% - 0.65% 

d2,*d2, = 

Nj 0... 0627 — 40.072 

N, = 

MWEFF 0.032% — 0.042% 

NalEFF = 0.152% 


kk 


# 


f is defined in equation IV.12. 


Ranges of values indicates that error varied 
angle, energy, or run. 


a RY far A 9 4 


with 


This total error is a worst case estimate obtained by 


adding the errors linearly. 


# 


100 


a 


7 — i ae 
Py — e _ - yp ; 
x >. a 5G “a 7 ‘ . ) 
— 2 = : a : 
>a Oe 7 = 
‘ fas : ‘ a a 
*esineutizganoy #3 hoe 4% al evos7a ye 
: . 7 ~~ 
some aan 
7 my 
2 rf g — 7 
= ait De 
Lr 7 
_ ae 
2 = 
Tr ‘ qv « ¥ 7 a ’ 
Salle \ 22 §3420- 
{ ‘c as =) a 
: = ‘ a 
5 : ie.¢ — Sates 
é é q 
= ae 


, ; r : 7 — y 
‘ a ; f 


~ : 
7 ri 4 
= 74% - Sab 
hs 
ae ha = 4.9 
aac 
‘vias P 
¥e0 = r4,.0 ° 
- co. = we-re’, Es 
4 a a * 
= 
. yn 


Pea = be. 
¥ 
ws 4 *e N 
2 é 5 v= Fic.9d - - 
e -, pa 
; a a . a A “et »? 
i . : 
4 7 ~~ 


; : 7 ’ a 

Sheva noianugs ab ie 

situ betaav Fetwe sed) ssset het Aleeki 
~fet es bs 


. 


boudesdo sdemijte a8A2 14308 2 “ % 
_ ¢ a ae oe watt vaot 


=. 
7 


CHAPTER V 


RLVLENG “OF THE) DATA 


The final stage of the analysis involved extracting 
unpolarized cross sections and analysing powers for the 
ONDE single particle valence states, (ld3,9 ,28 1/79, lds ;o); 
as a function of EDIF from the partially polarized cross 
sections stored in the ESUM versus EDIF histograms. The 
functional dependence of analysing power and unpolarized 
cross section on the partially polarized spin t and ¥ cross 
sections is given in equation AIII.10 of Appendix III. 

The computer software used for fitting were a TRIUMF 
program (OPDATA) and a CERN program (MINUITS). Notes on 
some specialized aspects of the fitting and on the 
calculation of the errors from the results of the fits are 
included in Appendix III. 

In a given 20 MeV EDIF interval, partially polarized 
cross sections for each of the valence states were obtained 
by fitting peaks to data over a portion of the ESUM range. 
Examples of ESUM spectra are shown in Figure V-l where the 
Epitmtesndeey datautorss0,:—-o4. 4for ‘channels’ 100.5 ‘to° 200.5 
in ESUM and EDIF = 20 MeV are presented. Since, as 
explained in Chapter IV, difficulties with the energy 
calibration resulted in small errors in summed energy which 
differed for each EDIF region, the amount of the shift from 
the correct value was determined in the fitting from the 


position of the ld3, 9 state. In Figure V-1, where the ESUM 
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EDIF=20 MeV. The positions of the oO Ge 
Valence™= states Cid 5775. 5287s de eee 
indicated along with those of possible Ca 
lp states. The wide and narrow-fit windows 
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scale is encotrrebkted for Stheriashift, eathee regions)! from 
channel 180.5 to 194.5 were the portions of these spectra 
Fit *inesthe final analysis. For tthis tdatival,es the t ofimad 
analysis determined the shift to be 0.3 MeV. Since the 
3.2 to 4.4 MeV ESUM resolution due to detector resolution 
and other experimental effects is large compared with the 
expected width of a particle stable state of less than an 
eiVi,¥ a Gaussian curve was taken to be an adequate 
representation of the peak shape. 

Spectroscopic factors derived from the analysis of a 
*0Ca(d, 2He) experiment (DO-76) suggest that most or all of 
ines astrengthi for) -eachs*of \tthe aaa and 28 1,,5 single 
particle states is concentrated in a single peak while that 
for the lds,» state is spread out over a number of peaks 
(over ~6 MeV in binding energy) as shown in Table VIII. 
The spectroscopic factors for the ld3,5, 281/2, and lds /o 
states’ (3.7, 1865, 4596) ‘are®"9323° -~83Z, and" ~83Z%0 of ‘their 
expected shell model values of (2J+1) respectively. The 
peak positions and spectroscopic factors from (DO-76) for 
these 3 states are fairly consistent with results of an 
earlier *Ga(p,d) experiment (MA-72) as well as with the 
results of a *Ca(e,e'p) experiment (MO-/76)- Consequently, 
our modelitors fitting stheadata utiddizged aesingle Gaussian 
peak for, eachseof thee ldes,, and 28,75 -S8tates- and a 
combination of 15 Gaussian peaks for the ld-,, state. For 
a given spin direction, there** were three’ variable 


amplitudes, one for each of the valence states (where the 
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amplitude for the ld,;,, was distributed on the basis of the 
spectroscopicy efactors in Table VIII). The relative 
positions of each of the peaks was fixed according to the 
excitation energies in Table VIII. 

The..shift in ESUM energy due to difficulties in the 
energy calibration was dealt with by allowing the positions 
of the peaks to shift up or down in ESUM with the spacing 
between the peaks remaining fixed. By varying the 
parameter for the position of the ld3,5 peak, XI, the 
fitting routine could determine the best position of the 
peak for the fit. The difference between the best value of 
XI in the fit and channel 191.67 was the assumed energy 
shift over the entire ESUM region fit. Spin t and ¥ data 
have identical energy calibration shifts and thus were fit 
together. 

The height for a single peak at a particular channel 
in ESUM, “XT; twas found by evaluating a function of the 
form 


(V.1) A*exp(-C52%*(XT-XI+EXC)7) , 


where A is the amplitude, 
XI is the position of the 1ld3,, peak, 
EXC is the excitation energy (relative to the 
ld3,5 peak) of the peak, 
C52 is a peak width parameter related to the 
standard deviation of the peak position, o, and the 


FWHM of the peak by 
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_— 
(V.2) C52 = = = i (FWHM in MeV) . 
(o} 


The area under a Gaussian is given by 


— Acyv¥t 
(V.3) Area = C52 : 


To remove the contribution of the 160(p,2p) lpj, 2 and lp3 92 
states from the *%Ca data, four peaks (two for each spin 
direction) were incorporated into the model as well. The 
values and statistical errors for the !°0 cross sections, 
fe one 
the results of a previous 200 MeV experiment (KI-80,KI-76). 


and analysing powers, nee were obtained from 


With the relations 
tye CO agri tPAl) 


do,( 160) a eRe) ato, + (oop ebay 2]h/? 


where POUSLOIM Spite, © SOP for spin ++, 


Spi ?0) 2 Aop( 1%) are the partially polarized 
cross sections and errors (ub/sr 2MeV) for. va 
polarization, P (where equation AIII.1l has been 
used), 
the 0 peaks “were, ancorpotrated into. the “fit. with 
amplitudes and errors given by 


(v.5) A(?&)+#4a( 180) = [o,+tAo, Jeg , 


be — 
. : ro 
7 
: —— ale 
s ? ps ‘ uae 2&2 v¢ yl “2 i. - e 
- J ‘S e -' ) = . °° = - * 
; - Ree 7 * = 
ety = 
7 
: Ren 
ie + 
J 


‘a sou ig @ ve isuaved a-se8ay ae 


‘ yeh “ 
° cz + 
7 2c2 = 
oo Hl el 
ia a VG aiid 30. A pr¢S@ tarao sid ~ 
4 jnsq ney sad 20° ans. mex? 
oe 


1132 oSat tsesegteqer  ste0 


7 : A : ; Pies: 
s 4 d4 ‘os e10T7e leads )isz2 hae pau 2 


- a oe. 
‘ ahha - eA & = “9 a a = %, 4 * hea * ‘ 
; fee A ee Hog. gpateyinns Sas ‘ 
eG sea ces Ms #% @OEF se0 tvs 36 ae $69 23 sab yt 
= ] 7 > a 8 
snetseten nts ast 
j * se pe pd Cc: «= cote: 
¥ . 
af ~ 
ps aK. 
| pany + Shae) Eg eeryh i= Cooke 
‘ ." . | re rn f ‘ ‘gs 


besJ og ¥ s 2 tag ee | ss iq Seea.** re : 
"oe 7 \ 
5o7 f "se \de) rset%9 Bas wataee asco) 


ssc est J.715A weltetpe 47 sk 2 3 ee 
<r owl a me a 7 


7 a Ld = 


date 317 eds oon) peaectaseank “snmih. sen amy 


107 


and with the overall normalization adjusted to the *lca 


datanby “ther facton 


C5226 RAGE 
(V-6) g =———— , 
YT 
where RALe esse ethem stra toe0 + 169 to 4+0ca atoms in the 
target. 


The angle pairs of data taken in this previous 
experiment did not completely overlap with all the angle 
pairs of the current experiment so that the fitting was 
performed only for energy differences between -100 and +100 
Mevecormehner-ane le pairs 040-94 , 47°-47°, 30°=30°, 30°35", 
Om A een 2 ee Oa Am SO 30° — 54° 5: and 30°-62°°. 
Results for three other sets of data where the 281/92 state 
was dominant (near conditions of zero recoil) were also 
obtained (67°-25° and 25°-67° (EDIF=120 MeV), 49°-35° and 
52°-40° (EDIF=40 MeV)). 

The amount of energy lost to recoil (~*0.3 MeV at 150 
MeV recoil for *Ca(p,2p) and ~0.4 MeV at 100 MeV/c recoil 
for 169(p,2p) (where the !°0 1p states have maximum cross 
section)) was small compared with the FWHM resolution of 
Bite CO, 4.4 MeV.. Thus, the shifting of the positions of 
peaks in ESUM due to varying recoil energy losses as the 
angles or energy differences changed was ignored in the 
model, and recoil energy losses were neglected in the 
positioning of the peaks. 


Even with the addition of these !°0 peaks, an '‘extra' 
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peak (B, for background) was required in the model to 
account for an additional feature in the spectrum at an 
excitation energy of 12.4 MeV. Without the inclusion of 
this peak, the x2 increased substantially. The nature of 
the peak that appears throughout the data at this energy is 
unknown® but one possibility jiswithat sitiuis eip,y75 strength 
from *O0ca(p,2p): This peak is part of the shoulder at 
approximacely eo channedwe)! 175)..5 4erarecls!] .Ose ine? theses pin~o14 
spectrum of Figure V-l. A second interesting feature 
appears between channels 160-5 and 178.5 (centered at 
FE69-5)— in@erhe.*spin + spectrum Pof* Fiscure “Vis This may 
possibly be lp3,,5 strength from *0Ca(p,2p). As this second 
feature is well outside the region of the data fitted, no 
account of it was taken in the model. Both of these 
features appear throughout the data in regions well away 
from 28,, state dominance and were most prominent in the 
30°-54° data which had the best statistics. 

The fitting of the data proceeded in two passes. In 
the first pass, the model was developed. A linear least 
squares routine (OPDATA) was used to minimize a x2 function 
of the form given in equation AIII.2. In summary, the 
model consisted of 40 peaks with a common width parameter. 
For each spin direction, there were seventeen 7° Ca peaks, 
two }&9 peaks, and a background peak. Ten parameters in 
the model were always free: four peak amplitudes for each 
spin direction (three for 40ca, one for the background 


peak), the energy calibration parameter (XI), and the 169 
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to «*%Ca sratio (RAT). The Fortran functions used to compute 
the value of the model for a particular spin direction are 
listed in Appendix IV where the 169 peaks are in a separate 
function from the rest. The same functions were used for 
both spin + and ¥ data. 

More *°Ca data (270%) were. acquired at 30°-54° 
(LF-RB) than for other angle pairs. For this reason, much 
of the early development of the model was performed using 
this data. As the background peak was most prominent in 
the 30°-54° data, it was used to determine the B peak 
position relative to the ld3,5 peak. The value obtained, 
12.4+0.3 MeV, was thereafter adopted as a fixed parameter 
in the model. 

As all the spin t+ and ¥ data for a given angle pair 
and detector combination should have the same resolution, 
the width parameter was established by first fitting 
several of the ESUM spectra from that data set with the 
width parameter, C52, free. By statistically averaging the 
result for C52 from three or four spectra according to 
equationgrAl bie354 ayevalvesr foreeC52 stand tit sceéeriror +t iwas 
established and all subsequent fits for that data set were 
performed with C52 fixed. These values and errors were 
also used in the second pass of the fitting. 

For a given spectrum, the range used in the fit was 
the same for the spin t+ and ¥ data and was chosen so as to 
include as much of the region from the B peak (on the low 
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possible without including other background features. The 
number of channels in these fitting regions (which we call 
the wide-fit window) ranged from 16 to 19 while the number 
of degrees of freedom for a fit varied from 22 to 28. The 
wide-fit limits (inclusive) for 30°-54°, EDIF=20 MeV are 
indicated by single-headed horizontal arrows in Figure V- 
ie 

All data were fit using the model and the wide-fit 
windows described above. The values of RAT, except for 
those from the EDIF=+1100 MeV spectra, where the statistics 
were generally quite poor, were used to determine the 169 
Pome. Caltration in thes 'iCa target. A histogram of these 
values of RAT is shown. in Figure V-2. This distribution 
resembles a broad Gaussian peak. However, the statistics 
are too poor to make any definite conclusion. A simple 
average and the sample standard deviation were used as our 
conservative estimate for RAT and its error. This yielded 
a value for 1°0/*%Ca of 0.32440.078 which was used in the 
second pass of the fitting. 

In the second pass (performed using the CERN fitting 
program, MINUITS (JA-75)), the model was basically the same 
as? s£0r ne first pass. tn weadgdiltion. sto othe ten -iree 
parameters of the first pass, C52 and the four 169 peak 
amplitudes were free parameters. However, RAT, C52, and 
the 166 amplitudes were constrained from deviating 
significantly from their previously determined average 


values. For each of these six parameters, a constraint 
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FIGURE V=2 
Distribution of RAT's, relative numbers of 
169 and *°%Ca atoms, produced from first pass 
imabyedschicye: “Weise Yale (eee c Thess centrodds fof sre 
distribution is indicated by an arrow. 
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term of the form given in equation AIII-.4 was added to the 
x2. In this manner, the uncertainties in these six 
parameters were reflected in final values for all the free 
parameters, their errors, and the x2. 

The upper channel limit for each second pass fit was 
Pheg same Saaatnat ifomethes corresponding first ‘pass. fit. 
However, the lower limit was chosen to be the channel 
mearesttoetll.> channels, below. the. position «of ,zthe dd, 75 
peak (XI) as determined from the first pass. In this 
fashion, the second pass ‘'narrow-fit' windows were chosen 
in a consistent manner for all of the data. In addition, 
while the size of this window ensured the inclusion of 
almost all of the lds,» state, it excluded a large fraction 
St the sbepeak resion thus minimizing the dependence of. the 
final results on the accuracy of the determination of the 
position of this peak. Varying the position of the B peak 
by amounts equal to its standard deviation (10.3 MeV) 
changed the. results. for,the, ld.» and .28 ;,5 states. by <24 
and had a negligible effect on the 1ld,,5 state. Also, the 
values for XI from the second pass fits differed by <0.5 
MeV from those of the first pass. The number of channels 
within these narrow-fit windows ranged from 14 to 17 and 
the number of degrees of freedom ranged from 19 to 25. The 
narrow-fit limits (inclusive) for 30°-54° EDIF=20 MeV are 
indicated by double headed horizontal arrows in Figure V- 
1. 


The result. of fitting the 30°-54° EDIF=20 MeV data is 
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shown in Figure V-3. The individual peaks are_ shown 
(solid lines for *0Cg peaks, dashed lines with circles for 
169 peaks, dashed line with +'s for background peak) along 
with the overall fit (dashed line) and the individual data 
points along with their errors (*'s with bars). Fits which 
produced negative amplitudes were redone with a penalty 
function added onto the yx* as shown in equation AIII.5. 
Five spectra gave slightly negative amplitudes C<S0ie> 
ub/sr*Mev) for the B or 23,5 peak and were ve-fit. 

From the *° a amplitudes, C52 values, and their errors 
as determined from the fits, the partially polarized cross 
sections and oe were calculated according to equation 
ALLL .6.< These values, in turn, were used to derive the 
unpolarized cross sections, analysing powers, and their 
errors by means of equations AIII.10 and AIII.11. The 
mesults for 29°-47°, 47°+29° (LF=RB, LB-RF), 30°-54°, 54°- 
SO aw sand), 25 —6/ 4, 67° —-25°. were averaged ‘statistically 
according to equation AIII.3. The final results for the 
o Ca unpolarized cross sections and analysing powers are 
presented and discussed in the next chapter. 

A general overview of the quality of the fits is 
illustrated in Figures V-4 and V-5. These figures contain 
the distributions of the percentage points and the we per 
degree of freedom respectively for all of the second pass 
bee Bey BP A percentage point is the fractional probability of 

2 


Obtaining, a qralue for -x- fequal,. to or jess than the»value 


obtained in the fit; its interpretation is discussed in 
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FIGURE V-3 


Results of fitting LE-RB=((30°-54—>). EDIF=20 
MeV data. Data points and their errors are 
indicated by asterisks and bars and the 
overall “fit “by aes dashed line. The *%Ga 
(143 79:7228 1/2, lde7o) andes “O8Glp oy tre 
peaks are shown as solid lines and dashed 
lines with circles respectively while the 
background peak is denoted by dashed lines 
with pluses (+). The energy calibration 
shift has been applied to the ESUM scale. 
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FIGURE V-4 é 
Distribution of percentage points from the 
second pass fitting of the data. 


FIGURE V-5 
Dist mi buttons soc xe per degree of freedom 
from the second pass fitting pot™ tnemdaea. 
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Appendix III. The percentage points are fairly evenly 
distributed with a small amount of skew toward the high 
end... (The x4 per degree of freedom histogram exhibits a x2- 
like distribution with a peak at a value of ~1.2. Neither 


distribution contains any disturbing features. 
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CHAPTER VI 


COMPARISON OF RESULTS WITH DWIA CALCULATIONS 


The values of the cross sections and analysing powers 
for the ld3; 2 (circles) and ld@’s;> (crosses) states obtained 
from the fitting are presented in Figure VI-l. For each 
angle pair, the data are plotted as a function of EDIF. In 
addition, the 2s8,,;,5 cross sections and analysing powers are 
plotted as a function of THETAp») Cthe>ip-p= scattering 
angle in the center of mass of the two protons) for 
conditions where the recoil momentum is at or near zero 
which includes: 25°-67°, EDIF=120 MeV (THETA))=52.1°); 
BUS 62° 5 EDIF=100° MeV “(62<1°); 30°-54°, EDIF=80 MeV 
C66... ):3 55a 49. TE DEEeAO, MeVea(s 7.7 -'): and. 940°=52.°, 
EDIF=40 MeV (81.5°). Away from zero recoil momentum, the 
2581/2 strength was small compared with that of the %=2 
states; consequently, its cross sections and analysing 
powers were poorly determined and are not shown. 

Eves Gurvesu in, Figure: Vil are the results of  DWIA 
calculations performed with the code THREEDEE (CH-79). 
This code performs the calculations according to equation. 
Abe and “employs va “direct S-dimensional numerical 
integration avoiding the use of angular momentum coupling 
coefficients. This represents a substantial improvement 
over previous codes which did not include the effects of 
spin-orbit distortion. The calculations were performed at 


the kinematic conditions for which the data were obtained 
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BLGUREM VL 1 

Cross sections and analysing powers for the 
ld375, ide7yo, “and 92s 75 states eResul Esme. 
the ld, 7os(cincles) “andm@id-) 53 (crosseq) mare 
presented as a function of energy sharing 
(EDIF) -for the angle (pairs “indicated mwhiwe 
the 2s 0/58 “nesultses at eee prec emt od aurd cummed 
function “of the, frees p—pmecenter scm mascs 
scattering angle (as calculated using a 
half-off-shell prescription). The DWIA 
calculations for the id= )5, 1d.) and else 
states are given as solid, dashed, and solid 
lines respectively while the dashed line in 
the Dillo © One 28,/2 analysing powers 
corresponds COs stn emer) Comm) meron qe ceen 
powers. The DWIA calculations were 
performed using spectroscopic factors one- 
half the simple shell model values otezs.). 
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and the dashed and solid curves are interpolations between 
the calculated values. Preetne ‘plots: fore the ld. 45" ‘and 
lds;yg states, solid and dashed curves are used respectively 
wove soltc™ curves sare;.used* in™*the plots tfor”’ the “2's;7>5 
state. ALI calculations “were *done*“with “an -optfeal 
potential including spin-orbit terms. Pn*Yvadditions>*the 
eased “tine vin thee pio © os ~- 2's") > Sanalysing**powers 
corresponds to the free p-p analysing powers resulting from 
excluding the spin-orbit terms in the optical potential. 
The cross sections are in units of ub/sr *MeV and, in order 
to plot the DWIA values and the measurements on the same 
scale, spectroscopic factors one half the shell model 
values of (2J+1) were used in the DWIA calculations. 

In the calculations, the bound state wave functions 
were calculated from a Woods-Saxon potential with spin- 
orbit terms. The radius, diffuseness, and well depth 
parameters are taken from values derived from fits to 
elastic electron scattering data as well as the binding 
energy of the particular single particle state (EL-67). 
ihe "binding energies*used in®the*’calculations for *the® Idiz\;5 
arte seq poe s cares aLhe= those =siven *in¥*lable ®VE white *a 
Varue "for 'the> ld state of 15.07 MeV (based on an average 
of the "binding energies "of ithe various ld), ‘components 
weighted by their spectroscopic factors) was used. As the 
calculations are not especially sensitive to small changes 
in the binding energy, this is an adequate approximation 


foresethe gtldic} es bind lig? “energy. The optical potential 
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(SC-80) used was derived from a global fit to angular 
distributions of p- +a elastic scattering cross sections 
and analysing powers including recent TRIUMF (200 MeV) and 
Indiana (80,135,160,181 MeV) data. This phenomenological 
optical potential utilized relativistic kinematics and a 
semi-relativistic Schrodinger-type equation. The potential 
does not contain an isospin-flip component which would 
result" Pfnvsa® chargettexchiange conmtrabutiion »tolcothie 
calculation. The functional form of this potential, its 
energy dependence, and the modifications necessary to 
implement it in the code are outlined in Appendix V. The 
difference between p- +a total reaction cross’ sections 
(o;) predicted using this optical potential and measured 
values is <9%Z at energies <140 MeV. Although the 
difference increases slightly to 14.5%  ati.181, MeV, 
distortion plays a relatively less important role in this 
energy region for an incident energy of 200 MeV. 

In the plane wave limit, the p-p interaction is half- 
off-shell. If an on-shell prescription is used, the half- 
off-shell nature of the interaction leaves some ambiguity 
as to the momentum at which to evaluate the p-p t-matrix. 
Two on-shell prescriptions are the initial and final state 
prescriptions. They correspond to evaluating the t-matrix 
at the momentum in the center of mass of the initial and 
gtruek Sprovons 2 Cinitialbivs tate) nor) <dnr’ that ..of. the “two 
outgoing protons (final state). More recently, a half-off- 


shell prescription in which the t-matrix is evaluated 
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between the initial and final state momentum (to account 
for the momentum and energy mismatch) has been incorporated 
into the calculations in a manner detailed in (MI-79). The 
tematrix Let factorized) intoiltwo® parts: the ‘on-shell t- 
matrix (in the form of partial wave amplitudes) which is 
evaluated from experimentally determined phase shifts (AR- 
77), and a real off-shell extension function calculated 
from the Mongan potential (M0-69) - a separable nonlocal 
phenomenological nucleon-nucleon potential. The results of 
calculations employing the initial and final state 
prescriptions differed slightly from those of the half-off- 
shell prescription by <+10% in cross section and <+0.08 in 
analysing power. All results shown correspond to the half- 
off-shell prescription. 

The modification of the free p-p interaction in the 
presence of the nuclear medium is manifested, in part, 
through off-shell behavior. theGelftect of “tthe “distorting 
potentials taking the interaction fully off the energy 
shell has been investigated (using the Mongan potential) 
for L69(2,2p) (MI-80) according to a method outlined in 
CRE- #3) Compared to the halftvotfi-sheli« calculations, 
there were significant changes only for forward angles 
(30°-30°) where the analysing powers are reduced by *<0..05 
and cross sections increased by <10%. 

petisecond effects of" the Cniucilear* meditin %on**PHe 
interaction is through Pauli blocking. Pauli prock ine 


refers to the suppression of amplitudes for p-p scattering 
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that leave one of the protons with a momentum equal to that 
offs Piomens off these protonerdof/) therlLcore. Calculations 
incorporating Pauli blocking performed by Miller and Thomas 
(MI-80) on !*0(~,2p) indicate a change of ~0.05 in 
analysing power and an increase of <l0Z in cross’ section 
compared with the half-off-shell calculations. Because 
they do not appear to affect the results of the 
calculations strongly, neither Pauli blocking or fully-off- 
shell calculations were performed in the present analysis. 
A careful comparison of the present DWIA calculations 
with the data reveals a number of trends. At equal proton 
angles and energies, the symmetry of the kinematics 
precludes + ithesgexistience\,, of -a,w preferned).direction. 
Consequently, the analysing power should and does vanish 
under such conditions in the data. Throughout the data, 
when the cross section calculations are normalized to the 
data, the resulting spectroscopic factors are ~50Z to 552% 
of the shell model values of (2J+1). In contrast, studies 
of other reactions (DO-/6,MA-72) have yielded spectroscopic 
factors lying in the range of the shell model values which 
is more consistent with naive expectations. Bo ren47 “tel, 
as the EDIF is varied, the recoil momentum reaches a 
minimum of ~60 MeV/c and here one expects a dip in the =2 
strength while the &=0 strength peaks. This is confirmed 
by both the data and the calculations. The agreement 
between the data and calculations for the &=2 analysing 


powers is quite good for forward angles (30°-30°, 30°-35°, 
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30°-40°) but deteriorates steadily as one or both angles 
increases in the backward direction. POPMecChe. (28.5 75 
analysing powers the difference between the DWIA 
predictions and the free p-p analysing powers is small and 
the data do not favor either set of values. The results of 
the previous !°0(~,2p) experiment at 200 MeV have been 
similarly compared with predictions from DWIA (KI-80). In 
this case the cross sections and analysing powers of the 
lpij,2 and 1p3,;5 states were calculated without spin-orbit 
distortion. The results were quite similar to the present 
case in that the cross sections were overestimated (by 
~80Z) and the agreement with the analysing powers was good 
POoCeeLOLrward wane les "but Spoor “as: one -or -both» angles 
increased. 

For 30°-30° and 30°-54° (£=2) and the 281, 7 cases, the 
calculations were redone with the real and imaginary spin- 
orbit terms in the optical potential set to zero. This 
caused an increase in the p-*%Ca Oy of <7% at the highest 
energy with no change at the lowest energy. The results 
are «shown ‘in Figure Vi-~-2. There were small but non- 
negligible changes in the (p,2p) cross sections and small 
changes in the analysing powers which produced no general 
improvement in the agreement with the measured analysing 
powers. 

The ov tical model potential describing the 
interactions when a proton passes through a nucleus 


can be calculated from the nucleon-nucleon (N-N) potential. 
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FIGURED Vi-2 
As in ftgure Viel except that the  sopin-ovba 


term has been turned Ontets in the DWIA 
calculations. 
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Pmt s Sicontext ,Téthe *antrinsicitnonlocality fof withe V N-N 
potential leads to a nonlocal optical potential. As well, 
nonlocality may be caused by antisymmetrization, as in the 
Hartree-Fock method of constructing self-consistent 
potentials for antisymmetric wave functions, or because the 
pEeatron? inh the elastic channel moves under the influence of 
coupling to the other channels (AU-65). 

Nonlocality is partly accounted for by means of an 
energy-dependent local optical potential, the calculated 
wave functions of which are asymptotically equivalent to 
those of its nonlocal counterpart. However, inside the 
nucleus, such local potential wave functions are not 
equivalent Bio the corresponding nonlocal potential 
wave functions. Calculations have shown that _ the 
amplitudes of wave functions derived from an attractive 
wonlocal.? potentialvrare isslenifiicantly » reduced; inside the 
nucleus compared with the corresponding amplitudes from a 
local potential which produces asymptotically equivalent 
wave functions. In the LPOG ¢2p) analysis, an 
investigation of this effect eas performed by multiplying 
each of the distorted wave function amplitudes by a factor 
dependent upon the nuclear matter density, a, where 
m is an effective mass characterizing the nonlocality 
of the potential. Values for this factor come from (JE-/77) 
as does the radial dependence of the nuclear density. This 


nromlocality  colmrection (method ds) was included in 


cabeullati ons ecfor a30°-30%9 B30 54°>e(R=2) Band, the 2847/2 
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values and the results are illustrated in Figure VI-3. 

A second method (ii) of adjusting the distorted wave 
functions to correct for nonlocality has been investigated 
using the WKB method (HO-80,AU-65,PE-63,PE-62) and is 
incorporated in the present DWIA code. In this method, the 
degree to which the amplitude of the wave function changes 
as a result of having a nonlocal component (such as the 
real central- part) in ‘the optical model potential is 
parametrized by a range parameter, 8. Wich>s 6 eset to 1.0 
fm and applying the correction corresponding to having a 
nonlocal real central potential to each of the distorted 
proton wave functions, calculations were performed for 30°- 
30°, 30°-54° and the 2s,,) cases. The results are shown in 
Figure VI-4. 

The two methods of incorporating nonlocality produce 
virtually indistinguishable analysing powers while the 
cross sections from method (i) are ~8%Z higher than those 
from method (ii). The nonlocal analysing powers differ by 
<0.03 from those of the local calculations. The nonlocal 
cross sections (method ii) are generally 20% to 50% lower 
than the corresponding local values. As B increases, the 
cross sections decrease and at a value of *1.15 fm the 
spectroscopic Ane for the =2 states, as obtained by 
normalizing the calculations to the data, increase to 
~(2J+1) in regions where the struck particle momentum is 
>55 MeV/c. In -the P°9(p,2p) analysis, method (i) 


produced effects similar to those observed here. 
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FIGURE VI-3 
As in figure VI-l except that a nonlocality 


correction (method i) has been incorporated 
in the DWDEA calculations. 
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FIGURE VI-4 
As in figure Viel excepts that a. nonloca ery 


correction (method ii) has been incorporated 
in the DWIA calculations. 
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While the inclusion of nonlocality in the calculations 
improves agreement with the cross sections, the inability 
to reproduce analysing powers for #0 states in certain 
kinematic regions remains. In response to these concerns, 
Miller (MI-80) has recently pointed out that the directions 
of the protons at the interaction point can be altered by 
Chenwjrefreacitviwve. citactspyotwethe opticalatpotenttiial . 
Consequently, the use of the asymptotic center of mass 
angle between the emerging protons in the calculation of 
the free p-p observables may be inappropriate. This, along 
with reaction localization in the target nucleus, may lead 
to a shift in analysing power the magnitude and direction 
of which would be correlated with kinematics and which 
would be the same for states of the same 2 ( 240) but 
different j. Preliminary studies of these effects indicate 
shifte insethe 1ds>5» and lds, 2 analysing powers producing 
closer agreement with the data. The effect on the 2817/2 
state is thought to be negligible. Detailed calculations 
are planned. 

For good shell model nuclei, if the spin-orbit term in 
the binding potential was zero so that the j= e41/2 and j=2- 
1/2 states were degenerate, then the net polarization of 
the subshell would be zero and the following relation 
between the effective polarizations, Porrf, of protons 
of the same & but different j would be expected to hold to 
a good approximation: 


CViis 17) PPI eda tak 20) Sel) Pe jms 2) : 
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The effective polarization is, in fact, insensitive to the 
details of the bound state wave function. Consequently, 
turning on the spin-orbit term in the binding potential 
should not affect the validity of equation VI.1 (MA-79a). 
In the absence of spin-orbit distortion, the effective 
polarization for a state may be obtained simply from its 
(p,2p) analysing power by means of equation AI.4. 
Using the free p-p analysing powers, A, and correlation 
parameters, CoN? obtained from phase shifts with a half- 
Vo euctieerrescription, othe. values of -RPogr for the 
tds75 and (jlde), states were determined. The resulting 
values of 2°P ogg (1d 3/2) (circles) and ~3°Po¢,¢(145/2) 
(crosses) are plotted as a function of energy sharing in 
Figure ViI-5. The results are consistent with equation 
VI.1 when both angles are in the forward direction. 
However, poorer agreement is observed as one or both angles 
increase in the backward direction. A similar analysis of 
the PS O(p 22'p) data (MA-79b) yielded comparable 


results. 
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CHAPTER VII 


CONCLUSIONS 


This experiment has yielded extensive information on 
the behavior of the cross sections and analysing powers of 
the ld3,5, 28 ,,/9, and ldc,» states of tigate The pattern of 
agreement between the DWIA calculations and the measured 
£=2 cross sections and analysing powers is remarkably 
similar to that observed in the earlier analysis of &=1 
states from the 169( 3, 2p) experiment. The inclusion of 
nonlocality in the DWIA calculations reduces the cross 
sections, allowing the spectroscopic factors to increase 
toward values expected of these nuclei. 

Although the inclusion of spin-orbit distortion in the 
present analysis has a marked effect upon the DWIA 
calculations, it does not significantly affect the quality 
of the agreement with the data. Bor ethei2sq ,5ostateyuthe 
data are consistent with either the inclusion or exclusion 
of spin-orbit dependence in the optical model potential. 

Evidence for the existence of 4%a lpi yo and 1p 3,5 
strength at ~21 MeV and ~30 MeV separation energy appears 
throughout the data although the energy resolution and 
statistical accuracy of our measurements in these regions 
is too poor to reach any definite conclusions. 

DWIA has proven capable of successfully predicting 
analysing powers for medium and light weight wees like 


*0¢a and 10 in restricted angular regions. ihe strong j- 
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dependencer'ofiethe analysing -power has been clearly 
demonstrated in these measurements. Given the current 
state of agreement between DWIA and experiment, (p,2p) 
reactions can be used as a spectroscopic tool to identify 
new pairs of states of the same & but different j provided 
one limits oneself to appropriate kinematic regions. The 
results of this experiment thus serve as an incentive for 
future efforts in both the theoretical and experimental 
investigation of quasi-elastic scattering with protons. 

A new measurement of tea (Beep) would help to 
resolve a number of the questions left outstanding in this 
experiment. With better ‘energy resolution (71 MeV), a 
quantitative investigation could be performed upon _ the 
feature which may be a "Ca IPipoeelaLe «aw DeLLer statistics 
should enhance the broad feature which may possibly be Stas 
lp3;2 strength. The analysing powers of these two features 
could then be compared with DWIA calculations to determine 
if they exhibit behavior consistent with lpiy2 and 1p3 79 
states. This identification would be facilitated by taking 
measurements at forward angles, for example 30°-30°, oe 
DWIA has demonstrated its ability to reproduce the data. 

The current experiment has stimulated the beginning of 
an examination of the effect of reaction localization 
coupled with refraction of the distorted waves in the 
mucleus on the calculations. Such investigations may 
extend the reliability of DWIA analysing powers to any 


kinematic region. 
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Comparisons of DWIA predictions with data are less 
“forgiving” when the data includes polarization data in 
that more stringent constraints are placed upon the theory. 
There is some indication that as the quality of the optical 
model information. inputs to e.the . DWIA, calculations, has 
improved, the implied spectroscopic factors have changed 
significantly (KI-80). Consequently, it might be wise to 
ne-examing, earlier results, in .this light. Only when DWIA 
is able to predict analysing powers and cross’ sections 
simultaneously will one be confident in the resulting 


nuclear structure information. 
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APPENDIX .I 


DWIA DESCRIPTION OF (,2p) 


The. description sof jquasi-—elastic (p,2p), reactions with 
the Distorted Wave Impulse Approximation (DWIA) has been 
reviewed by Jacob and Maris (JA-66,JA-73a). Their notation 
is shown, in, Figure,.I-2. Assumptions made in DWIA 
treatments of quasi-elastic scattering include the 
following: 

i) there is a single violent interaction between the 

incident proton and struck nucleon and there are no 

other strong interactions between the incident proton 
and the residual nucleus, 

5 | the violent interaction can be related to free 

proton-proton scattering, 

pe A OB the distortion factors (defined as the ratio of 

plane wave functions to distorted wave functions) do 

nots chance japprmectably. overpethe;-urange,. lof, this 
interaction, 

( ii) and iii) are collectively a statement of the 

impulse approximation) 

iv) the core interacts with the incoming and outgoing 

nucleons as a refractive and absorptive medium. This 

ieeoescra bed ine Lerner of, an optical model. potential 
which distorts the wave functions of the incoming and 
outgoing protons, 


v) the single particle shell model wave functions are 
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taken as an adequate representation of the overlap 

integral between the target and residual nuclei. 

In the present context of generating distorted waves 
for protons with energies from 50 MeV to 200 MeV, the best 
optical model potentials available are empirically obtained 
from fitting elastic cross section and polarization data. 
A more desirable approach would be to generate optical 
potentials from multiple scattering theory, such as KMT 
CKE=39)°3 utilizing the fundamental nucleon-nucleon 
interaction. However, at the current time, such an 
approach has not yielded optical potentials for energies 
less than *400 MeV which are consistent with the elastic 
cross section and polarization data (NA-80). 

The differential cross section under the assumption of 
no spin-orbit dependence in the optical potential was shown 
by Jacob and Maris to be given by a product of a factor (K) 
involving the phase space and spectroscopic factors, the 
distorted momentum distribution of the struck nucleon (G), 


and the free proton-proton scattering cross section, 
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d3o es PRS 
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lf there is a single value for the initial orbital angular 
momentum quantum number, &, of the struck proton, then G is 


given by 
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The Wy are weight factors involving Clebsch Gordan 

coefficients. With ‘the further assumption® of zero 

distortion (the plane wave limit), eo nee eet eae our ler 

transform of the initial bound state wave function of the 
m 

struck proton, 9: 


nm ae > n> > 
(AI.1.2) Gy = fexp@iieg) iechien, (cdr ; 


> 
Here Koy is the momentum of the recoil nucleus. 
In spherical coordinates, the momentum dependence of 
Gh enters viap’a spherical ‘Bessel)" function, jg(k,_47) 


(ME-66). For small values of the arguement of jg, 


Q 


(ATS. 3 ) A-1%) . 


di peee (Cle 
Consequently, at zero recoil momentum, the momentum 
distribution for %#0 states goes to zero while that for 2=0 
states is at a maximum since jg is at a maximun. Ea 
practice, the minimum in the momentum distribution for #0 
states at zero recoil is somewhat filled in due to the 
effects of distortion and the finite angular acceptance of 
the detectors. 

The cross section for scattering free protons with 


parallel or antiparallel spins is given by 


dort dort 
es a, cae e + e e e 
(AT.2) d 2 d 2% {2 CES ee? ae ware Spe ‘ 


Here P and P fe are the polarizations of the incident and 
struck protons respectively while A and Cnn are the 


free EG p-p tanalysing~-ppowers#® and «correlation *tarameters 
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respectively. Prtetet fie tand eA aredemeaisumed with 
respect to an axis perpendicular to the scattering plane 
and the polarizations are >0 for spin +t, <0 for spin. +. 
Z2Lquintermediigteggenergics, aCynolislalarge dand 
positive. Thus the cross section for protons with parallel 
spins is 3 to 5 times larger than that for protons with 
antiparallel spins. Bound protons with the same & (#0) 
but different j can be effectively polarized in opposite 
directions in quasi-elastic reactions. Expression (AI.1) 
leads one to expect states of the same & but different j to 
exnibiee different (behavior in their (p,2p) analysing 
powers. — For %£=0 states, the effective polarization 
vanishes, and thus the (f,2p) analysing power reduces 
to that for free p-p scattering. 
The (p,2p) analysing power, AysitOt a) particular 


state is obtained from the spin + and spin + cross sections 


by 
dv s(t) d3ca (+) 
did Q5dE did 25dE 
GAL) 2Aa = : 
decom (xh ie comet) 
d 21d 25dE d2)dQ5dE 


Substituting equations AI.1l and AI.2 into AI.3 and solving 


for Perf, one obtains 


aie — A 
(AI.4) P oft = ie. SS She ° 
NN y 


Thus, by measuring the analysing power, Ay» for a state 


and given the free p-p values A and Cun? the effective 
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polarization of the struck nucleon can be determined. 

The inclusion of spin-orbit dependence in the optical 
model potential seems to be important in fitting elastic 
and reaction cross section as well as polarization data 
(NA-80). Lf included, Sepin-orbit distortion ‘upsets. the 
factorization of the cross section and expression (AI.1) is 
no longer valid. Rather, as shown by Jackson (JA-76), the 
cross, aeection. is ‘factorized in the amplitudes and 


calculations become considerably more complicated, 


AI.5) Rea |r 
Sea df, dQ5dE ine if 


Here C is a phase space factor and Sij vepresents the 
appropriate sum and average over initial (i) and final (f) 


states. In DWIA, Try is given by 


yey 1 
6 = De e Z e e 
Ga. 63, Wa. tino, § 5s 3,6 J) (J,M,IM|I,M,) ( tmzo, | JM) 


* Ca -—* eal Ss = (3 + Ca es 
q or eX: t q es te r aX: ! q ,ar r 
Wo 2 Wy Hy 1 RJ 55% 0 
> 
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1 + Le oe 
aa ae ee ae (pie 5 
139'K5 Ky) [tp loo os pCR gtk, 1)? , 


¥ 


where kek oks ak yey are the lab momenta of the incident 


proton, the outgoing protons and the recoil nucleus, 


> > > 
dg241,2d5 are the momenta in the Atl center of mass 


(JASG5)—_ conjucate sto, theyespatial coordinates 
appearing in the Schrodinger equation for the 
distorted wave functions, 

1 > > 


76k, 7k) are the relative momenta of the 


a eee 
76 Ko ee 


protons in their cem.- system before and after 
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scattering, 
(...-|--) are Clebsch Gordan coefficients, 


ey y6*9) are fractional parentage coefficients, 
if. 


Oo,u are spin quantum numbers, 
sete are the distorted wave functions, 
N is the spectroscopic factor corresponding to 
the final state (=2J+l in the single particle shell 
model where J (or j) is the total angular momentum 
quantum number of the state), 
and a=(A-1)/A. 
Omitting isospin, the overlap integral ae of the 


initial and final nuclear states is given by 


M > 
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Here s,; and 6, are the spin quantum numbers for the struck 
proton and Xe is a spinor. Livectnewecontext. Sof “oir 
*[Ca(p,2p) calculations, the overlap integral is taken to 
be a pure single particle state with a single value of 2& 


ANG od. Thus, all configurational mixing is ignored and 


there is only a single non-zero value of ey ae 
ae 
AD Swipe St auement wabout .the result of the 
complication arising from the inclusion of the spin-orbit 
term is difficult to make, but one can make a few 


generalizations. The relative ‘sizes “of various p-p 
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amplitudes determine the polarization of the final state 
protons resulting from the free elastic scattering. This 
polarization of the particles after the interaction can 
affect the probability of their emergence from the nucleus 
resulting in a selection ‘of specific p-p. amplitudes. 
Hence, the -apin-orbit distortion results in a modification 
Grew loesceu relative atrengsths “through ‘the differing 
absorption and distortion of particles with opposite 
polarizations. injvecdadt tion, the polarization of ‘the 


incoming proton may be affected strongly (SC-78). 
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APPENDIX II 


SPECTRA STORED IN MEMORY DURING EXPERIMENT 


Number Parameter Length 
(Words) 

1 Y-LF 128 
2 X-LF 128 
o Y-RF 128 
4 X-RF 128 
i X-LB 128 
6 Y-—LB 126 
7 X-RB 128 
8 Y-RB 128 
9 TDC 256 
10 TDC 256 
Lt TDC 256 
132 EDG 256 
te ADC : 256 
14 ADC 256 
PS ADC 256 
16 ADC 256 
a | EUG 256 
18 TDC 256 
19 ADC + ADC a2 
20 ADC + ADC le 
21 ADC VS ADC 4096 
22 ADC VS ADC 4096 
Pace ADC VS sADC = 1024 
24 ADC VS ADC 1024 
26 DCR LO 16 
26 DCR HI 16 
Pas Scaler 768 
28 ADC 256 
29 ADC pasa s 


Ch) S-Pevnoe puleer 


Gating Comment 
Requirement 


LRP 

LF °=-P 

RF °-P 

RF *-P 

DR =P 

PoP 

RB*-P 

RB *-P 

(LF*RB)*-P LF Start,Stop=RB 
(LB*RF)*-P LB Start,Stop=RF 


2 Left Start,Stop=rf 
ae Right Start,Stop=rf 
LF NaI-LF 

RF NaI-RF 

LB NaI-LB 

RB NaI-RB 


(LF *RF) e-P RF Start,Stop=LF 
(LB *RB) *-P RB Start,Stop=LB 


(LF RB) NaI-LF + NaI-RB 
(LB *RF) NaI-LB + NaI-RF 
(LF °RB) NaI-LF vs NaI-RB 
(LB °RF) NaI-RF vs Nal-LB 
(LF ¢RF) NaI-LF vs NaI-RF 
(LB *RB) NaI-LB vs NaI-RB 
Time 
LF Plastic=-LF 
RF Plastic=RF 
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APPENDIX III 
STATISTICAL TREATMENT OF ERRORS, STATISTICAL AVERAGING, 


DEFINITION AND INTERPRETATION OF x* 


Counting statistics in a spectrum are given by 
Poisson statistics which become Gaussian when the number of 
counts geper? bin, ce becomes’ large (e220) - The 
condition, e220, Wasensatistied”.fors alle dataleanalysed 
in gthis,,experimenty-and hence all statistical quantities 
were assumed to follow Gaussian statistics. Yc was taken 
as an adequate estimate of the statistical error in the 
measured number of counts, .c. 

Consider a function of n parameters, F(P1>Poo***oPL): 
Some of the parameters will have corresponding 
uncertainties and may be correlated with other parameters. 
Still others will be pre-determined and have no uncertainty 
and no correlation with other parameters. The uncertainty 


in F is given by 


n n n 
oF 2 oF oF 
XSL IIE pe lk AF = ys eer 1 ay aie ee aN 
( ) eee a aan ers ig ieee Pert ney ee, 
ak ob i 
oie er nnd er terms [1/2 
where the summations extend over the parameters which are 


not pre-determined, 


Ap, is the error in the ith parameter, 


A(p.p.) is the covariance between the ith and jth 
ij 


parameters. 


In general, a mathematical model which is to be fit 
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to a .set of measurements, y ay, (i=l,m), must make a 


i 
preciction, Fy 


» corresponding to each measurement, yjy- 
Thespredictions come ‘from, the evaluation of a function, 
FCP) >Pos*++sPL)> Some of the parameters will be determined 
and thus fixed while others are free to be varied by a 
fitting program. The purpose of the fitting program is to 


find a set of values for the free parameters which minimize 


the value of a function called x2. In the simplest case, 


x- may be computed from 
m Fa Coee pes... Sprj=y, 
tie , pte A ges n tey2 
CATTIC2Y x 4 ae eee : 
i 
2. 


This expression for yx“, when applied to a spectrum 
where m is the number of bins, is approximately true if the 
number of counts is not too small (>5). Strictly speaking, 
the yy should be independent, normally distributed random 
variables and Ay; must be their standard deviations. 

As a special case, if a series of measurements, 
ee (i=1,m), are made of a single quantity, y, the 


statistically weighted average for y and the corresponding 


error are given by 


m 
x Po hy Os = 
Ye WAES Dre —= a és } 
Peis -2 
epee’ pg? 
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In addition to returning a set of values for the 
parameters which minimize x2, the fitting routine also 
returns 4 symmetric metrix variously referred to as the 
eovariance, variance.) or) error matrix. The ith diagonal 
element ot this mateix=is the variance (or the square of 
Cucmeseebaatds deviartony mot, the (ith “free parameter. 
Changing one of the parameters by an amount equal to its 
standard deviation and refitting with this parameter fixed 
would result in x2 increasing by 1. These standard 
deviations are taken as the errors in the parameters, 
Api. The (i,j)th element of the matrix is the covariance 
between the ith and jth parameters, ACPiPy)- 

Pneaihe Scourge ote developing.. a .fittine even a 


parameter may be used which has previously been determined 


independently to a certain accuracy, Pi td: In this case, 


it may be desirable to let the parameter remain free while 
constraining it with this additional information. This 
2 


conetraintwmisweimplemented siby vincluding ia® térm cin the Lx 


Sune tion, 


rela CLES he 


Another situation that arises is the desire to restrict the 
value of a parameter to a particular region. This is done 
by means of adding a penalty term into the x2. An example 
of a x* function with both a constraint term and a penalty 


term is 
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Here the third term (the penalty term) helps to constrain 
the dth parameter to non-negative values. Penalty terms 
mayo (helpeito, prevent. the fit from producing ‘unphysical 
values for parameters such as negative amplitudes for 
peaks. 

The value of yx? may be used to determine the degree 
to which the data was successfully fit by the model. One 
measure for this is the value of x2 per degree of freedom. 
The number of degrees of freedom is given by the number of 
data points minus the number of unconstrained free 
parameters. A second (not independent) measure based on 
the x? and number of degrees of freedom is called the 
percentage point. Pablesttotf* percentage ‘joints are 
available in standard mathematical handbooks (BE-76). Its 
value ranges from 0 to 1 and it represents the fractional 
probability of obtaining a value for x ¢ equal to or less 
than the value of x? observed. Fou, ayrsing lel sfit. <a 
definite conclusion based on either of these two quantities 
is difficult to outline but very general statements can be 
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made. Fits with  a>>¥ per degree of freedom (percentage 


HOM.) Vor yl (70525 indicate a reasonably jg00d fit; 
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wery® Parngeivalues » >48(70.295), bindicate atpoor fit and put 
the validity of the model into question; very small 
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too good and May esuegsest i fhatarthe errors in the data 
have been overestimated. 

For a number<of sets. of data fit by a single model, 
if the model describes the data and the data and variances 
follow a normal distribution the histogram of percentage 
points should be flat. 

The model used in the fitting involved Gaussian peaks 
StmecrewmtOruerciven in vequatton Vv. 1. The result of fitting a 
set of data were two amplitudes for each peak in the model 
(one for each spin direction). In addition, for some fits, 
tie width, C52, was.a variable parameter. Let 
At+AAt be the amplitude for the spin +t peak, 

AYtAAY be the amplitude for the spin ¥ peak, 

C52+AC52 be the width parameter, 

1A spat og be siie beam polarizations of the spin -+,% .data 
(where both values are between O and 1). 

In terms of these parameters, the partially polarized spin 


t and *¥ cross sections and their errors are given by 
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where equations V.3 and AIII.1 have been used. 
Let of and of represent the cross sections for 100Z 
polarizations. The relations between the partially 


polarized and 100% cross sections are 


CALET . 2) es 5 °(14P +) +04 + 5°(1-P +) +04 
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< 
i] 


1 
5*(1-P 4) sot Peec(ltPt) sot |. 


Solving for ot, oO we have 


(1+P +) eot - (1-Pt) cot 
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1) ame Sy 


IFFT) Ot — Ol =P) sot 
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Pt++P + 


The unpolarized cross section and analysing powers are 


defined as 


(AIII.9) oe THe ; 
l= ot-—oy 
y ~ ot+or J 


and in terms of the partially polarized results, these 
quantities become 
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Pt tatey, ron e€quationa ALLL. «ALLI.10 and recalling that 


oe and om are independent quantities, the errors in Oo, and 


Ay are given by 


ee : 2 ; 2 11/2 

(AIII.11) 409 = oaaay [(P + dot) + (Pt dot) ] 
iW SE athe [CoyeAot) 2 4+ (otsAos) 21/2, 

y Lie sei 


2 
(P¥*eot+P toot) 
P P 


The statistical errors in Pt and P¥ are sufficiently small 
to be ignored. The scale error due to the polarimeter p-p 


calibration does not affect the fits. 
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APPENDIX IV 


FITTING FUNCTIONS 


THIS CALCULATES THE CONTRIBUTION OF 40CA(P,2P) TO THE 
ENERGY SUM SPECTRUM AT THE POINT XT. 

A52,A32,AS12,AB ARE THE AMPLITUDES OF THE 40CA 
1D52,1D32,2S12, AND BACKGROUND PEAKS. 

SH IS A PARAMETER USED TO SEPARATE SPIN UP AND DOWN 
DATA. 


+ ee FF HEH He HF HH HF HF KH HK He F 


FUNCTION FITCAS2 SA32,AS12),AB ;C52,BB,X1,XT,SH) 

IMPLICIT REAL*8 (A-H,0-Z) 

DATA DSUM/4.96/ 

CSQ=-C52*C52 

CONST=XT-SH-XI 

FIT=A52*0.91/DSUM*DEXP(CSQ*(CONST+5.27)**2) 
+A52*0.64/DSUM*DEXP(CSQ*(CONST+5.61)**2) 
+A52*1.25/DSUM*DEXP(CSQ*(CONST+6.34)**2) 
+A52*0.10/DSUM*DEXP(CSQ*( CONST+6.50)**2) 
+A52*0.10/DSUM*DEXP(CSQ*(CONST+6.77)**2) 
+A52*0.18/DSUM*DEXP(CSQ*(CONST+6.96)**2) 
+A52*0.10/DSUM*DEXP(CSQ*(CONST+7 .20)**2) 
+A52*0.47/DSUM*DEXP(CSQ*(CONST+7 .43)**2) 
+A52*0.10/DSUM*DEXP(CSQ*(CONST+7.78)**2) 
+A52%*0.35/DSUM*DEXP(CSQ*(CONST+8.17)**2) 
+A52*0.24/DSUM* DEXP(CSQ*(CONST+8.43)**2) 
+A52*0.22/DSUM*DEXP(CSQ*(CONST+8.55)**2) 
+A52*0.10/DSUM*DEXP(CSQ*(CONST+8 .90)**2) 
+A52*0.10/DSUM*DEXP(CSQ*(CONST+9.10)*¥*2) 
+A52*0.10/DSUM*DEXP(CSQ*( CONST+9.75)**2) 
+A32*DEXP(CSQ*( CONST) **2) 
+AS12*DEXP(CSQ*(CONST+2.52)**2) 
+AB*DEXP(CSQ*(CONST+12.4)**2) 

RETURN 

END 


THIS CALCULATES THE CONTRIBUTION OF 160(P,2P) TO THE 
ENERGY SUM SPECTRUM AT THE POINT XT. 

APE2Z, AP32-ARE THE AMPLITUDES OF THE 160. 1P12° AND 1P32 
PEAKS. 


* 


FUNCTION CTM(AP12,AP32,C52,XI,XT, SH) 

IMPLICIT REAL*8 (A-H,O-Z) 

CSQ=-C52*C52 

CONST=XT-SH-XI 

CTM=AP12*DEXP(CSQ*(CONST+3.797)**2) 
+AP32*DEXP(CSQ*(CONST+10.11)**2) 

RETURN 

END 
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APPENDIX V 


OPTICAL MODEL POTENTIAL FOR p-‘*°Ca SYSTEM 


The optical potential used in the calculations 
consisted of a Coulomb term, a complex valued nuclear 
eentral oterms, .and-~a complex valued nuclear Spin-or pit 
term 


CAVeL) AU Cr) =2.U - VE o(r,R5,a 9) 


coulis. 


ilw. - 4a Ww, ee + 
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So so oO) So 


> > 
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ie lie 
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where the form factors, f, have Woods-Saxon shapes, 


Be 
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ees R)/a 


The values of the various parameters as a function of 
energy are listed in Table X. The details of the analysis 
method used to derive this potential are given in (NA-80). 
The analysis used to obtain the optical potential was 
based on a semi-relativistic Schrodinger-type wave equation 
derived from the Dirac equation for a fermion moving in a 
Locair.z2eq scentral potential, U(r). As the spin-orbit term 
arising from the Dirac equation is negligible, the one used 
in this analysis is purely phenomenological. With the 
projectile Streated trelativistically and «the target. noa- 
relativistically, the radial wave equation for the relative 


motion of the system is given by 
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TABLE X 


Optical Model Values* for p- *%Ca System 
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mass of the p- *%a system. 


The corresponding non-relativistic Schrodinger 


equation used in the DWIA program is 
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Thus, in order to use the potentials in a manner consistent 
with the way they were derived, the following substitutions 


were made for the potential, V, and the reduced mass, u, 
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